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1.  INTRODUCTION 

The  use  of  microprocessors  in  advanced  control  schemes  has  resulted  in  an 
overall  improvement  in  the  engineer's  ability  to  monitor  and  process  information. 
Utilizing  the  microprocessor  requires  interfaces  between  the  system  to  be  con¬ 
trolled  and  the  controlling  processor.  If  the  plant  to  be  controlled  is  mechanical, 
mechanical -electrical  transducers  are  required  for  feedback  information  and 
electromechanical  transducers  are  required  for  output  actuation.  This  report 
deals  with  the  introduction  of  a  new  electromechanical  transducer.  The  medium  of 
mechanical  actuation  is  air;  therefore,  the  device  discussed  herein  is  an 
electropneumatic  converter. 

Pneumatic  systems  offer  the  advantages  of  high  power,  fast  response,  and 
economical  control  hardware.  The  electropneumatic  converters  that  have  been  de¬ 
veloped  so  far  have  tended  to  be  quite  slow  and  expensive.  What  is  needed  is  a 
low  cost,  low  electrical  power  consumption,  high  speed  transducer  to  take 
advantage  of  the  benefits  of  low  cost  and  reliability  which  could  be  obtained 
with  an  electropneumatic  control  system. 

An  electropneumatic  system  can  be  a  low  cost  control  system.  Because  of 
advances  in  semiconductor  technology,  electrical  signals  can  now  be  processed 
extensively  at  a  low  cost.  The  technology  to  convert  air  under  pressure  to 
mechanical  actuation  is  a  refined  art,  and  the  hardware  is  available  and  inexpensive. 
Therefore,  a  requirement  for  the  electropneumatic  transducer  is  to  devise  a 
minimunncost  device  to  keep  the  total  system  cost  low. 

Another  requirement  for  the  converter  is  reliability.  Reliability  is  an 
inherentjcharacteristic  of  pneumatic  actuating  hardware.  Therefore  reliability 
of  the  converter  is  required  to  avoid  the  possibility  of  a  weak  link  in  the  system. 
This  can  be  achieved  through  simplicity  of  design  and  a  minimum  of  moving  parts. 

Low  electrical  power  consumption  is  a  desirable  feature  for  a  new  design. 

If  the  power  consumption  is  low  enough,  the  converter  could  be  attached  directly 
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to  the  output  of  an  inexpensive  low-power  electrical  signal  source.  This 
electrical  signal  source  could  be  the  digital  to  analog  output  port  of  a  micro¬ 
processor.  This  would  eliminate  the  need  for  expensive  electrical  power  ampli¬ 
fication.  An  additional  benefit  of  low  electrical  power  consumption  is  an 
application  where  an  onboard  power  supply  is  required.  A  device  which  uses  low 
power  could  save  additional  weight  and  expense  through  the  use  of  a  small  power 
supply. 

A  final  design  consideration  is  speed  of  response.  The  microprocessor  used 
to  control  the  system  is  a  very  high-speed  device.  Pneumatic  to  mechanical  hardware 
is  capable  of  fast  response  when  the  volumes  are  small  and  air  pressure  is  high. 

The  interface  device  should  also  be  as  fast  as  possible  to  maximize  the  control 
capabilities  of  the  processor. 

The  general  design  plan  for  a  system  that  fulfills  these  requirements  is  now 
given.  A  low  electrical  power  consumption  electromechanical  component  is  required. 
The  mechanical  output  of  this  component  will  then  be  converted  to  a  low-level 
pressure  signal.  This  low-level  pressure  signal  will  then  be  amplified  to  a 
usable  level . 

The  device  chosen  for  the  electromechanical  conversion  is  a  piezoelectric 

bender  element  of  a  bimorph  construction^.  The  mechanical -pneumatic  conversion 

uses  the  bimorph  as  a  flapper  centered  between  two  nozzles.  The  device  chosen 

2 

for  the  pneumatic  ampl i fication  is  the  laminar  proportional  amplifier  . 

A  piezoelectric  bimorph  provides  a  displacement  for  a  given  applied  voltage. 
This  steady-state  relationship  between  the  voltage  input  to  the  bimorph  and  its 
output  displacement  is  nearly  proportional .  The  bimorph  has  very  low  electrical 
power  consumption  and  a  fast  response.  The  reliability  of  the  bimorph  appears 

^Application  Note,  Piezoelectric  Bender  Elements,  Piezo  Products  Division, 

Gulton  Industries,  Inc.  (1978). 

2 

F.M.  Manion  and  T.M.  Drzewiecki,  Analytical  Design  of  Laminar  Proportional 
Amplifiers,  Froc.  HDL  Sta te-of-the-Art  Symposium,  Vol.  1,  Adelphi,  MD  (Oct.  1974). 
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to  be  quite  good.  And  lastly,  in  quantity,  bimorphs  are  relatively  inexpensive. 
This  general  information  is  encouraging  enough  to  explore  incorporating  a  piezo¬ 
electric  bimorph  as  the  major  moving-part,  electromechanical  component  in  an 
electropneumatic  converter. 

While  the  piezoelectric  bimorph  can  output  a  steady-state  displacement  for  an 
applied  voltage,  its  output  force  capability  is  limited.  Static  and  dynamic 
fluid  flow  forces,  due  to  the  air  flowing  in  the  nozzles,  should  be  minimized. 

This  minimizes  the  load  forces  which  would  interfere  with  the  displacement  due  to 
the  applied  «oltage.  As  a  result,  a  low  supply  pressure  and  small  nozzle  diameters 
are  necessary  to  minimize  flow  forces.  This  results  in  small  pressures  produced 
downstream  of  the  nozzles.  These  low-level  pressure  differences  can  be  amplified 
to  a  usable  level  using  laminar  proportional  amplifiers. 

The  laminar  proportional  amplifier  (LPA)  is  a  no-moving-part  pneumatic 
amplifier.  It  is  capable  of  high  gain  and  high  bandwidth.  The  LPA  can  also 
operate  a  very  low  supply  pressures.  Because  it  operates  in  the  laminar  region, 
the  LPA  has  high  output  fidelity  and  low  noise.  The  use  of  a  piezoelectric 
bimorph,  functioning  as  a  flapper  between  two  nozzles  in  conjunction  with  cascaded 
LPA  amplifiers,  to  amplify  the  low  level  pressure  difference  available  down¬ 
stream  of  the  nozzles,  results  in  an  attractive  solution  to  the  electropi.eumatic 
conversion  problem. 

2.  CONCEPT  OF  VALVE  DESIGN 

3  4 

Figure  1  is  a  schematic  drawing  of  the  proposed  valve  concepts  ’  .  The 
valve  uses  a  piezoelectric  bimorph  element  centered  between  two  nozzles.  This 

^C.K.  Taft  and  B.M.  Herrick,  A  Proportional  Electro-Fluidic  Pneumatic  Valve 
Design,  20th  Anniversary  of  Fluidics  Symposium,  The  American  Society  of 
Mechanical  Engineers,  Chicago,  Illinois  (1980). 

4 

C.K.  Taft  and  B.M.  Herrick,  A  Proportional  Piezoelectric  Electro-Fluid  Pneumatic 
Servovalve  Design,  1981  Joint  Automatic  Control  Conference,  The  American  Society 
of  Mechanical  Engineers,  Charlottesville,  Virginia  (1981). 


Figure  1.  Flapper  Valve  Schematic 


flapper-valve  arrangement  is  pressurized  in  a  chamber  by  a  supply  pressure  P^. 

The  pressures,  P|  and  Pg,  downstream  of  the  nozzles  depend  on  the  location  of  the 
bimorph  relative  to  the  nozzles. 

With  a  voltage  applied  to  the  bimorph  as  shown  In  Figure  2,  there  is  a  result¬ 
ing  piezoelectric  effect  which  causes  the  bimorph  to  deform  in  such  a  manner  that 
it  will  be  closer  to  nozzle  2.  The  new  position  oft the  bimorph  results  in  a 
difference  between  the  restriction  at  the  entrance  to  each  of  the  two  nozzles. 

The  deflection  shown  in  Figure  2  causes  the  flow  area  at  the  entrance  to  nozzle  1 
to  increase  and  the  flow  area  at  the  entrance  to  nozzle  2  to  decrease.  This  will 
cause  the  pressure  downstream  of  nozzle  1  to  be  higher  than  the  pressure  downstream 
of  nozzle  2.  Therefore,  there  Is  a  low-level  pneumatic  signal  available  at  the 
output  of  the  nozzles.  The  pneumatic  pressure  difference  is  proportional  to  the 
electrical  input  signal  for  biirorph  to  nozzle  distances  less  than  about  one  eighth 
of  the  nozzle  diameter. 

This  pressure  difference  is  then  applied  across  the  input  ports  of  an  LPA 
While  the  flapper-nozzle  configuration  produces  a  small  pressure  difference.  Mis 
difference  can  be  amplified  by  using  several  of  the  LPA  amplifiers  cascaded  together. 

Figure  3  illustrates  how  an  LPA  operates^ — Ajet-of  fluid  flows  from  the 
supply  nozzles  and  traverses  a  distance  where  it  is  separate  by  a  splitting  wedge. 
With  zero  pressure  difference  applied  to  the  control  ports,  the  jet  is  undeflected 
and  assumes  the  centered  position.  The  same  amount  of  fluid  enters  each  output 
port  which  then  results  in  a  zero  output  pressure  difference.  When  a  small 
pressure  difference  is  applied  across  the  jet  in  the  vicinity  of  the  control 
region,  the  jet  will  deflect.  The  deflection  will  cause  a  pressure  difference  at 
the  output  ports.  The  gain  of  the  amplifier  is  controlled  by  its  internal  geometry. 
That  is,  a  small  jet  deflection  at  the  supply  nozzle  can  become  substantial  at 
some  distance  downstream.  Typical  state-of-the-art  amplifier  designs  result  in  a 
pressure  gain  of  ten  when  the  output  ports  are  blocked.  These  amplifiers  can 
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then  be  cascaded  in  series  and  parallel  for  Impressive  system  gains  in  pressure 

2 

and  flow  . 

The  piezoelectric  bimorph  and  the  LPA's  must  now  be  combined  in  some  way  to 
meet  a  set  of  design  requirements.  For  a  given  application  there  are  several 
factors  that  can  be  varied  to  change  the  valve  system's  output  characteristics  to 
suit  the  application.  Examples  of  output  characteristics  which  can  be  varied 
are  system  sensitivity,  bandwidth,  pressure,  and  flow.  These  factors  are  all  a 
function  of  the  size  and  physical  characteristics  of  the  piezoelectric  bimorph 
chosen  and  the  size,  element  thickness,  and  supply  pressure  used  in  LPA  assembly. 

The  LPA  is  made  up  of  a  series  of  laminations  stacked  on  top  of  each  other. 
The  number  of  laminations  used  in  each  stage  of  amplification  defines  an 
important  design  parameter,  the  supply  nozzle  aspect  ratio.  The  aspect  ratio  is 
the  supply-nozzle  height  divided  by  the  supply-nozzle  width.  Varying  the  aspect 
ratio  will  change  the  static  and  dynamic  response  of  the  amplifier  stage.  The 
aspect  ratio  is  one  of  the  easiest  design  parameters  to  vary  and  which  defines 
the  amplifier's  characteristics. 

Certain  qualitative  design  objectives  for  an  electropneumatic  converter  can 
be  noted.  A  given  application  will  require  certain  output  pressure  and  flow 
levels.  These  output  characteristics  are  a  function  of  the  LPA  design.  A  matched 
set  of  complimentary  LPA  stages  is  called  an  LPA  gain  block.  The  gain  block  may 
be  assembled  with  each  stage  arranged  in  series  with  the  stage  before  it.  It 
is  also  possible  to  arrange  a  gain  block  design  with  some  of  the  stages  in 
parallel  with  each  other. 

The  series  design  will  maximize  the  pressure  gain  that  can  be  obtained. 
However,  the  flow  level  will  be  relatively  low.  This  will  be  due  to  the  use  of 
amplifier  stages  with  successively  smaller  aspect  ratios.  The  use  of  a  smaller 

2 

F.M.  Manion  and  T.M.  Orzewiecki,  Analytical  Design  of  Laminar  Proportional 
Amplifiers,  Proc.  HDL  State-of-the-Art  Symposium,  Vol.  1,  Adelphi,  MD  (Oct.  1974). 
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aspect  ratio  at  each  level  of  amplification  is  desirable  for  proper  ampliier 
input-output  impedance  matching  and  to  maintain  laminar  flow  with  in  the  amplifier 
stages  that  operate  at  higher  supply  pressures. 

A  series  design,  with  some  or  all  of  the  amplifier  stages  within  the  gain 
block  consisting  of  identical  amplifiers  in  parallel,  will  increase  the  flow 
through  the  system.  This  particular  staging  technique  is  shown  schematically  in 
Figure  4.  The  complete  set  of  amplifiers  comprising  stage  1,  stage  2,  and  stage  3 
are  all  in  series  with  each  other.  Amplifiers  Ug  and  X£  are  in  parallel  with 
each  other,  and  the  amplifiers  U^,  X^,  Yj  and  Z2  are  also  in  parallel  with  each 
other.  The  system  of  Figure  4  will  deliver  a  higher  flow  at  any  given  output 
pressure  when  compared  to  a  similar  system  which  is  arranged  without  the  additional 
parallel  amplifiers;  but  this  will  occur  at  the  expense  of  pressure  gain.  This 
reduction  in  pressure  gain  is  due  to  decreased  amplifier  input  resistance  between 
each  series  stage. 

In  addition  to  designing  an  overall  gain  block  which  is  capable  of  delivering 
the  required  output  flow  and  pressure  levels,  each  amplifier  in  each  stage  must 
be  designed  separately  to  match  the  amplifier  before  and  after  it.  This  is 
accomplished  by  impedance  matching  each  stage  of  amplification  along  with  matching 
jet  deflection  from  one  stage  to  the  next.  This  matching  is  necessary  to  assure 
that  in  all  stages  the  laminar  jet  will  sweep  the  same  angle  for  a  given  input. 

This  will  prevent  premature  saturation  of  the  gain  block  system. 

Also,  each  stage  of  amplification  must  be  checked  for  speed  of  response. 
Considering  the  requirements  discussed  above,  one  would  usually  desire  to  build  a 
system  which  is  as  fast  as  possible.  This  is  accomplished  through  proper  choice  of 
both  the  IPA  gain  block  and  the  piezoelectric  bimorph.  A  smaller  bimorph  generally 
results  in  a  faster  transduction  of  electrical  energy  to  motion  at  the  end  of  the 
bimorph.  As  the  length  of  a  bimorph  is  reduced,  the  motion  of  its  end  diminishes, 
and  overall  system  sensitivity  is  reduced.  This  may  be  compensated  for  by 


’•  v 


additional  stages  of  amplification,  but  will  result  in  reduced  overall  amplifier 
gain  block  bandwidth.  This  reduction  in  bandwidth  is  caused  by  increased  signal 
path  length  and  lower  bandwidth  found  in  low  supply  pressure  amplifiers  which  are 
required  to  amplify  the  reduced  output  at  the  nozzles.  There  is  an  optimum 
obtained  by  trading  off  bimorph  size  for  the  number  of  amplifier  stages. 

Finally  the  static  and  dynamic  fluid  forces  have  to  be  calculated.  Any 
significant  force  on  the  end  of  the  bimorph  will  degrade  its  motion.  A  maximum 
dynamic  flow  force  coupled  with  a  maximum  static  pressure  area  force  has  to  be 
determined.  Thus  based  on  a  tolerable  maximum  flow  force,  a  corresponding  maximum 
supply  pressure  to  the  bimorph  can  be  selected.  The  bimorph  supply  pressure  must 
then  be  less  than  this  maximum  value.  j 

The  above  design  requirements  demonstrate  the  need  for  a  comprehensive  mathe¬ 
matical  model  describing  the  behavior  of  the  valve  system.  With  a  sufficiently 
descriptive  model,  a  designer  could  take! a  set  of  specific  design  requirements 
and  use  the  model  to  decide  if  a  given  b|imorph-LPA  system  can  meet  those  require¬ 
ments.  This  approach  to  problem  solving  would  both  save  time  and  give  the  designer 

j 

more  insight  into  how  the  system  works. 

The  total  valve  system  is  broken  down  into  two  subsystems.  The  first 
subsystem  is  the  electromechanical  system  consisting  of  the  electrical  connection 
to  the  bimorph  and  the  piezoelectric  bimorph  itself.  The  second  subsystem  is  the 
mechanical -pneumatic  system,  which  consists  of  the  nozzles  and  the  LPA  gain  block. 
The  coupling  between  the  two  systems  is  the  motion  of  the  end  of  the  bimorph  and 
the  fluid  forces  on  it.  Section  3  will  deal  with  the  development  and  verification 
of  the  electromechanical  model.  Section  4  will  similarly  deal  with  the  development 
of  the  fluidic  model.  These  two  models  can  then  be  used  with  computer  simulation 
to  develop  a  prototype  design  which  will  then  be  built  and  tested. 
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3.  DEVELOPMENT  OF  THE  ELECTROMECHANICAL  MODEL 

A  dynamic  model  describing  the  output  motion  of  the  end  of  the  piezoelectric 
bimorph  for  a  given  electrical  input  is  required.  This  model  will  enable  the 
designer  to  select  a  bimorph  and  determine  if  it  is  suitable  for  the  design 
requirements  that  are  presented.  Before  a  model  can  be  developed,  some  general 
background  information  on  piezoelectric  materials  and  their  properties  is  necessary.. 
This  will  make  the  different  components  of  the  model  easier  to  understand. 

3.1  Piezoelectric  Background  Material 

Certain  materials  generate  an  electrical  charge  when  they  are  deformed.- 
This  effect  was  first  discovered  by  the  Curie  brothers  in  1880.  It  was  also 
discovered  that  this  is  a  reversible  effect  and,  therefore,  if  a  charge  is  applied 
to  one  of  these  materials,  deformation  will  occur.  Basic  research  in  this  field 
was  performed  by  Cady5,  who  determined  that  the  electrical  polarization  of  the 
substance  was  proportional  to  its  strain  and  that  the  effect  was  sign  sensitive. 

This  effect  was  called  the  piezoelectric  effect,  and  the  materials  which  exhibited 
this  effect  were  called  piezoelectric  materials.  It  was  determined  that  the  piezo¬ 
electric  effect  could  only  occur  in  materials  that  had  anisotropic  crystal 
structures. 

The  piezoelectric  effect  must  not  be  confused  with  a  similar  but  fundamentally 
different  effect  called  electrostriction.  The  electrostrictlon  effect  occurs  In 
all  solid  dielectric  materials.  It  is  a  reversible  effect,  which  makes  It 
similar  to  the  piezoelectric  effect.  The  difference  occurs  because  the  electrical 
polarization  caused  by  electrostriction  is  proportional  to  the  square  of  the  strain,. 
This  causes  the  effect  to  be  sign  insensitive.  Electrostriction  is  also  smaller 
in  magnitude  than  piezoelectricity5. 

5W.G.  Cady,  Piezoelectricity,  McGraw  Hill  Book  Co.,  New  York,  1946. 

5S.Y.  Lee,  Piezoelectric  Actuators  for  Fluid  Control  Applications,  Eng.  Proc. , 

Fluid  Control  Systems,  Pennsylvania  State  Univ.,  (July  1965),  p.  45. 


The  two  groups  of  piezoelectric  materials  which  are  the  most  significant  and 
useful  are  the  natural  and  synthetic  crystal  group  and  the  polarized  ferroelectric 
ceramic  group.  The  crystal  group  is  comprised  of  materials  such  as  quartz,  lithium 
sulfate,  and  rochelle  salt.  The  ferroelectric  ceramic  gorup  is  comprised  of 
materials  such  as  quartz,  lithium  sulfate,  and  rochelle  salt.  The  ferroelectric 
ceramic  gorup  is  comprised  of  materials  such  as  barium  titanate  and  lead  zirconate 
titanate.  Before  any  comparisons  can  be  made  of  any  specific  examples  within  these 
groups,  definitions  of  the  general  piezoelectric  parameters  are  necessary. 

Two  main  families  of  constants  pertinent  to  the  piezoelectric  effect  are 
the  g  coefficients  and  the  d  coefficients.  The  g  coefficient  is  the  ratio  of  the 
electrical  field  strength  produced  divided  by  the  stress  applied.  The  d  co¬ 
efficient  is  the  ratio  of  the  electrical  charge  generated  divided  by  the  applied 
force.  The  nomenclature  developed  for  these  coefficients  requires  a  two-digit 
subscript.  The  convention  used  states  that  the  first  subscript  should  indicate 
the  electrical  signal  field  direction  and  the  second  subscript  the  stress 
direction.  Convention  also  states  that  the  3-direction  is  parallel  to  the 
direction  of  piezoelectric  polarization  of  the  material.  Figure  5  illustrates 
the  coordinate  directions  and  rotations7.  Since  normal  and  shear  stresses  can  be 
applied,  six  directions  are  specified. 

The  constant  that  relates  the  d  coefficient  to  the  g  coefficient  is  called 
the  absolute  dielectric  constant,  e.  The  absolute  dielectric  constant  is  the 
ratio  of  the  g  coefficient  divided  by  the  corresponding  d  coefficient.  The 
magnitude  of  the  dielectric  constant  gives  an  indication  of  the  electrical 
capacitance  and  dielectric  properties  of  the  material. 

The  product  of  the  g  coefficient,  the  material's  elastic  modulus,  and  the  d 

o 

coefficient  is  the  dimensionless  parameter,  K  ,  which  is  defined  as  the  square  of 
the  coupling  coefficient.  The  coupling  coefficient  is  a  measure  of  the  electrical 

7E.G.  Doeblin,  Measurement  Systems,  McGraw  Hill  Book  Co.,  New  York  (1975). 
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energy  available  for  a  mechanical  input.  The  magnitude  of  the  coupling  coefficient 
is  indicative  of  the  relative  efficiency  of  various  piezoelectric  materials. 

The  piezoelectric  effect  is  relatively  independent  of  temperature.  There  is 
a  temperature,  though,  above  which  the  piezoelectric  effect  disappears.  This 
temperature  is  called  the  Curie  point.  The  Curie  point  generally  occurs  at  a 
temperature  which  results  in  a  crystalline  change  of  phase.  This  breaks  down  the 
anisotropic  crystalline  structure  observed  in  piezoelectric  materials.  It  is  this 
change  of  phase  which  destroys  the  material’s  ability  to  exhibit  the  piezoelectric 
effect®. 

Now  that  general  piezoelectric  parameters  have  been  defined,  several  materials 
may  be  examined  for  their  desirability  as  the  material  which  will  comprise  the 
electromechanical  component  of  the  converter.  Table  1  lists  several  piezoelectric 
materials  and  their  properties. 

TABLE  1.  PIEZOELECTRIC  PARAMETERS  FOR  SELECTED  MATERIALS8 


Coefficient 

xlO'12 

meters 

volt 


Coefficient 

xlO"3 

volt-meter 

newton 


Dielectric 

Constant 

xlO*11 

farad 

meter 


E  K 

Elastic  Coupling 

Modulus  Coefficient 

xlO10 

newton 

meter 


Curie 

Pt. 

°C 


Rochelle 

Salt 

-165 

-93 

200 

1.93 

0.5| 

1  45 

Quartz 

-2.3 

-5.8 

4.5 

8.0 

o.  ic; 

1  550 

I 

Lithium 

Sulpha te 

-16 

-175 

10.3 

4.6 

0.36 

» 

Bari  urn 

Titana te 

-148 

-16 

1150 

11.1 

0.51 

)  125 

Lead  Zirconate 
Titanate 

-180 

-11 

2000 

6.0 

0.34 

330 

Q 

W.P.  Mason,  Piezoelectric  Crystals  and  Their  Application  to  Ultrasonics, 
Van  Nostrand,  New  York,  (1950). 
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Examining  the  d  coefficient  provides  an  opportunity  to  compare  the  motion 
output  for  an  applied  voltage.  Rochelle  salt,  barium  titanate,  and  lead  zirconate 
titanate  all  have  relatively  high  d  coefficient  values,  indicating  they  would 
provide  a  large  displacement  for  a  given  applied  voltage.  This  property  would 
be  advantageous  in  a  flapper  design. 

The  dielectric  constant  for  barium  titanate  and  lead  zieconate  titanate  is 
relatively  high.  Since  the  dielectric  constant  is  a  measure  of  a  material’s 
resistance  of  rupturing  with  an  applied  electrical  field,  both  barium  titanate 
and  lead  zieconate  titanate  would  have  relatively  high  allowable  input  voltages 
which  would  result  in  larger  displacements. 

Finally  the  maximum  operating  temperatures  of  the  materials  are  examined.  A 
high  maximum  operating  temperature  is  desirable  for  many  applications.  Normal 
military  specifications,  for  example,  require  a  maximum  operating  temperature  of 
83C .  Quartz  shows  a  maximum  operating  temperature  of  550C.  This  is  more  than 
sufficient  for  most  applications.  Lead  zirconate  titanate  has  a  Curie  point  of 
330C,  which  would  make  it  also  acceptable  in  most  applications.  Barium  titanate 
has  a  Curie  point  of  125C,  which  would  make  it  a  borderline  material  in  military 
appl ications.  Rochelle  salt  with  a  Curie  point  of  45C  makes  it  unacceptable  for 
military  applications. 

Considering  their  sensitivity,  dielectric  constant  and  Curie  point,  barium 
titanate,  lead  zirconate  titanate,  and  similar  materials  are  chosen  by  piezoelectric 
manufacturers  as  the  materials  for  production  of  piezeoelectric  manufacturers  as 
the  materials  for  production  of  piezoelectric  motors.  Barium  titanate  and  lead 
zirconate  titanate  are  members  of  the  polarized  ferroelectric  ceramic  group  of 
piezoelectric  materials.  This  group  is  different  than  the  crystal  group  in  that  a 
processing  procedure  must  be  used  to  give  them  their  piezoelectric  properties. 

The  word  ferroelectric  is  derived  from  a  dielectric  analog  Y  with  ferromagnetic. 
This  will  become  obvious  after  the  processing  procedure  of  the  ferroelectric 
ceramic  group  is  explained.  Ferroelectric  ceramic  substances  are  polarized  by 
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applying  a  strong  electrical  field  while  heating  the  substance  above  its  Curie 

point  and  then  cooling  it  below  the  Curie  point  with  the  electrical  field  still 
q 

applied  .  The  Curie  point  in  this  particular  case  is  the  temperature  at  which 
a  polarized  ceramic  loses  its  uniform  polarization  and  assumes  random  polarization. 
At  this  temperature  there  is  a  breakdown  in  the  piezoelectric  properties  of  the 
material. 

It  is  not  necessary  to  raise  the  temperature  of  a  ferroelectric  ceramic  to 
its  Curie  point  to  polarize  it.  However,  fully  polarizing  a  ferroelectric  ceramic 
requires  a  much  lower  electrical  field  strength  if  the  temperature  of  the  material 
is  raised  above  its  Curie  point.  The  polarization  process  is  more  easily  understood 
if  the  ceramic  is  considered  to  be  comprised  of  a  set  of  randomly  polarized  domains. 
Therefore  an  electric  field  applied  to  the  ceramic  will  result  in  an  alignment  of 
these  domains.  A  stronger  electric  field  will  cause  more  alignment  to  occur. 

There  is  also  a  strain  and  resulting  change  in  dimension  that  occurs  as  a  result 
of  the  applied  field.  The  observed  strain  of  a  nonpolarized  ferroelectric  ceramic 
is  proportional  to  the  square  of  the  applied  electric  field.  Since  the  strain  is 
also  sign  insensitive  to  the  applied  electric  field,  the  effect  observed  here  is 
electrostriction  and  not  the  piezoelectric  effect.  Since  the  ceramic  retains  some 
polarization  after  an  electric  field  has  been  applied  to  it,  some  hysteresis  in  the 
steady-state  relationship  of  electrical  input  versus  displacement  output  should  be 
observed.  This  is  analogous  with  any  type  of  hysteresis  observed  in  ferromagnetic 
materials  and  the  reason  the  material  is  called  ferroelectric^. 

Once  an  electrostrictive  material  is  polarized,  it  develops  a  sensitivity  to 
bidirectional  signals.  Also  the  observed  strain  is  proportional  to  the  applied 
field  in  each  single  domain  of  the  ferroelectric  ceramic  crystal.  This  then 

9  ‘ 

E.O.  Doeblin,  System  Dynamics:  Modeling  and  Response,  Charles  E.  Merrill 
Publishing  Co.,  Columbus,  OH  (1972). 

^S.Y.  Lee,  Piezoelectric  Actuators  for  Fluid  Control  Applications,  Eng.  Proc. , 

Fluid  Control  Systems,  Pennsylvania  State  Univ.,  (July  1965),  p.  45. 
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indicates  that  a  pure  piezoelectric  effect  is  observed  in  these  polarized  materials 
This  results  in  input-output  characteristics  similar  to  those  found  in  more 
traditional  piezoelectric  materials.  Because  of  this  similarity,  the  same  modeling 
equations  that  govern  the  characteristics  of  piezoelectric  materials  can  be 
applied  to  ferroelectric  ceramics. 

It  is  now  desired  to  produce  a  maximum  displacement  for  a  given  applied 
voltage.  This  can  be  accomplished  if  the  bimorph  method  of  construction  is  used. 
Figure  6  is  a  cross-section  of  a  piezoelectric  material  with  a  bimorph  construction. 
Two  sheets  of  material  which  exhibit  piezoelectric  properties  are  each  plated  on 
one  side  with  a  layer  of  nickel  or  silver.  The  sheets  are  then  bonded  to  a  brass 
shim! to  form  a  bimorph.  This  same  concept  may  also  be  extended  to  multilayered 
piezoelectric  sheets  resulting  in  a  multimorph.  The  multimorph  will  produce  a 
larger  displacement  for  an  applied  voltage  but  consumes  more  power. 

Jhe  opposing  polarities  of  the  two  sheets  of  piezoelectric  material,  indicated 
I  ■ 

by  Figure  6  demonstrates  that  this  bimorph  is  to  be  connected  in  series  with  the 

i 

poweri  supply.  If  the  polarities  of  the  two  sheets  were  aligned,  then  the  bimorph 
would  be  connected  in  parallel  to  the  power  supply.  A  series  connection  requires 

i 

twice;  the  voltage  to  output  the  same  force  or  displacement  as  a  parallel  connection 
but  only  half  the  charge.  For  this  reason  a  bimorph  connected  in  series  to  the 
power  supply  will  have  a  steady-state  input  impedance  of  four  times  the  impedance 
of  the  parallel  connection. 

To  further  increase  the  output  motion,  the  bimorph  is  configured  as  a  bender 
with  the  piezoelectric  expansion  in  the  length  expander  mode.  The  pertinent 
piezoelectric  coefficients  in  this  mode  of  operation  are  the  g^  and  the  d^  co¬ 
efficients.  Therefore,  when  a  voltage  is  applied  to  the  bender,  one  piezoelectric 
component  expands  while  the  other  contracts.  This  causes  the  composite  to  bend 


r.‘  .  Anderson,  Theory  Ferroelectric  Behavior  of  Barium  Titanate,  Ceramic  Age, 
57',;,  (1951),  pp.  29-30,  33,  53-55. 


like  a  bimetal ic  strip.  This  bender  is  then  used  as  a  flapper  in  a  flapper-nozzle 
design. 

An  electromechanical  model  of  the  bimorph  bender  can  now  be  developed  and 
analyzed.  The  bimorph  can  be  represented  by  an  equivalent  mechanical  circuit  and 
an  equivalent  electrical  circuit.  The  two  circuits  are  then  combined  to  provide 
an  overall  model  relating  the  motion  of  the  end  of  the  bimorph  to  the  electrical 
input. 

3.2  Equivalent  Mechanical  Circuit 

The  first  mode  transverse  motion  of  the  bimorph  is  modeled  as  an  ideal 
spring-mass-damper  system.  Figure  7  shows  the  bimorph's  equivalent  mechanica' 
circuit  assuming  that  the  motion  of  the  end  of  the  bimorph  is  nearly  in  a  straight 
line.  This  model  has  been  found  to  be  sufficiently  accurate  for  the  first  im.de  of 
vibration  of  the  bimorph. 

Summing  the  forces  on  the  mass  element  in  the  X-direction  gives 

MX(t)  =  -K$X ( t)  -  BX  (t)  -  C-, V(t)  +  Fx(t)  .  (1) 

The  mass,  M,  of  the  system  is  an  equivalent  mass  found  by  matching  the 
frequency  of  the  first  mode  of  vibration  for  a  cantilever  mounted  beam  to  that  of 
the  simple  system  shown  in  Figure  7.  The  spring  constant,  K$,  is  a  parameter 
which  relates  the  force  necessary  to  produce  a  given  displacement  at  the  end  of 
the  bimorph.  The  damping  constant,  B,  describes  a  force  proportional  to  the  velocity 
of  the  end  of  the  bimorph.  The  product,  C^V,  is  a  coupling  force  between  the 
mechanical  circuit  and  the  electrical  circuit.  Finally  ( t )  is  an  external  forcing 
function  to  represent  the  fluid  flow  forces  acting  on  the  bimorph  in  the  X-direction. 

Taking  the  Laplace  transform  of  equation  (1),  setting  the  initial  conditions 


equal  to  zero,  and  solving  for  X(S)  yields 
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3.3  Equivalent  Electrical  Circuit 

Figure  8  shows  an  equivalent  electrical  circuit  for  the  bimorph  driven 
by  a  voltage  source.  The  bimorph  is  modeled  as  a  capacitor,  Cy,  with  a  leakage 
resistance,  Ry,  between  the  electrical  contacts.  The  current  source,  K^,  pro¬ 
portional  to  the  velocity  of  the  end  of  the  bimorph,  is  a  result  of  the  charge 
produced  when  a  piezoelectric  device  is  deformed.  The  product  is  the  coupling 
between  the  electrical  system  and  the  mechanical  system.  It  is  assumed  that  the 
system  is  driven  by  an  external  voltage  source,  VA,  along  with  an  external  series 
resistor,  Rft.  The  series  resistor  will  be  used  to  modify  the  frequency  response 
characteristics  of  the  bimorph. 

Summing  the  currents  at  node  V  gives 

K  x  (t)  -  c  imi  .  mi  .  . : . 1^1 * 0  oi 

v  m  I  dt  Ry  Ra  u  * 

Taking  the  Laplace  transform  of  equation  (3)  setting  the  initial  conditions 
equal  to  zero,  and  so’ving  for  V(S)  gives 


V  ( S )  = 


K  SX(S)  +  -~ 


Va(5) 


cts  +  \r +  r 

1  ,  rt  ra 


Substituting  equation  (4)  into  equation  (2)  and  solving  for  X(S)  results  in 


ctrs  +  T 

Fy(S)  -  -  VJs) 


-yV"Rr II 


,  CtMR  ,  CtDR  +  M  +  C,KjR  +  B 

s3  ,T_  +s2  .T_  r.  +  $ 


where 


ra  *  RI 


(6) 


3.4  Experimental  Verification 

Experimental  verification  of  the  model  was  desired  to  determine  its 
accuracy  and  resulting  suitability  for  design  work.  A  photo  interrupter  was  used 
to  measure  the  deflection  of  the  end  of  the  bimorph  resulting  from  an  electrical 
input.  The  bimorph  used  to  verify  the  model  was  of  a  larger  size  than  the  element 
which  was  eventually  used  in  the  valve.  It  was  necessary  to  determine  the  electro 
mechanical  system  parameters,  Cp  M,  K$,  B,  Cp  Kq,  to  properly  apply  the  model. 
The  physical  and  material  properties  of  the  test  bimorph  are  given  in  Table  2. 
These  properties  were  used  to  predict  the  electromechanical  system  parameters. 

TABLE  2.  PHYSICAL  AND  MATERIAL  PROPERTIES  OF  LEAD  ZIRCONATE  TITANATE  CERAMIC 

BIMORPH  USED  FOR  MODEL  TESTING 


Property 
Length,  L 
Width,  W 
Thickness,  T 

Nickel  plate  thickness,  tj 
Lead  zirconate  titanate  thickness,  t ^ 
Brass  shim  thickness,  t-j 
Density,  t> 

Piezoelectric  charge 
coefficient,  d^ 

Piezoelectric  voltage 
coefficient,  g.^ 

Nickel  elastic  modulus,  E, 

(see  Appendix  A) 

Lead  zirconate  titanate 
elastic  modulus,  E£ 

Brass  elastic  modulus,  E^ 


Measured  Value 
4.3  x  10‘2  m 
1.27  x  10"2  m 
4.9  x  10-4  m 
1.52  x  10'6  m 
1.93  x  10‘4  m 
1 .03  x  10”4  m 
5.8  x  103  kg/m3 

-1.8  x  10'10  m/v 

-1.1  x  10“2  v-m/N 

1.26  x  1011  N/m2 

1.09  x  1011  N/m2 
9.0  x  1010  N/m2 


■A'.'.J'v'  .W  ..nv J.s1.- ■  * '  W.V/ •'. 1 V. 
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Lead  zirconate  ti tana te  was  chosen  as  the  ferroelectric  ceramic  material 
used  to  comprise  the  bimorph.  It  was  chosen  because  of  its  high  d.^  coefficient, 
high  dielectric  constant,  moderately  high  coupling  coefficient,  high  maximum 
operating  temperature,  and  commercial  availability  in  bimorph  bender  form.  A  series 
connection  was  chosen  because  it  allowed  a  simple  electrical  connection. 

The  value  of  the  system  spring  constant.  K$,  can  be  predicted  using  beam 
theory  from  an  expression  relating  the  deflection,  a,  at  the  free  end  of  a 
cantilever  mounted  beam  to  a  force,  F',  applied  at  its  end. 


(7) 


This  can  be  used  to  obtain  an  analytical  expression  for  the  system  spring 
constant. 


(8) 


Figure  6  illustrates  that  the  material  comprising  the  bimorph  is  not  uniform. 
Instead,  there  are  laminations  of  material  running  lengthwise  down  the  bimorph. 

An  equivalent  system  spring  constant  is  calculated  by  creating  an  equivalent  beam 
with  a  uniform  modulus  of  elasticity  and  a  resulting  equivalent  width.  Equation  (9) 
can  then  be  applied  to  obtain  the  spring  constant. 

The  equivalent  width  of  a  lamination  varies  proportionally  with  the  material's 
elastic  modulus.  Therefore, 


Ei 

wE1  =  W^)  , 

Ep 

WE2 =  w(r}  •  and 

WE3  *  W^E'^  ’ 


(9) 

(10) 
(ID 


The  moment  of  area  of  the  beam  of  Figure  6  is  given  by 
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TI  WE3*3  +  2  3^2*2  +  WE2t2  2  + 


2  3  WE1t]  +  WE1t1[  2  +  t2Jj 


Substituting  the  properties  of  Table  2  into  equations  (6), (9), (10), (11),  and  (12) 

results  in  an  analytically  determined  system  spring  constant  of  *  310  N/m. 

The  system  spring  constant  can  also  be  measured  to  verify  this  calculation. 

Recalling  equation  (1),  the  spring  constant,  K$,  can  be  determined  by  measuring 

the  displacement  due  to  a  steady-state  force  at  the  end  of  the  bimorph  with  the 

,  •* 

terminals  short  circuited.  The  steady-state  conditions  would  require  that  X  *  X 
=  0,  and  short  circuiting  the  terminals  would  make  V  =  0;  therefore,  equation  (1) 
becomes  K$X  - .  This  can  be  measured  experimentally  using  a  photo  interrupter 
to  measure  displacement  and  a  Celsco  force  transducer  to  measure  force.  This  re¬ 
sulted  in  an  experimentally  determined  system  spring  constant  of,  K$  *  290  N/m. 

Another  system  parameter,  the  effective  mass  of  the  bimorph  is  obtained  by 
solving  the  partial  differential  equation  for  the  first  resonant  frequency  of  a 
cantilever  mounted  beam  modeled  by  a  distributed  parameter  technique.  This  first 
mode  natural  frequency  is  set  equal  to  the  natural  frequency  obtained  with  the 
single  lumped  parameter  approximation.  The  bimorph's  effective  mass  can  then  be 
found  in  terms  of  known  system  parameters. 

The  first  mode  natural  frequency  of  a  cantilever  mounted  beam  modeled  by  a 
distributed  parameter  technique  is  given  by  Oen  Hartog^  as 


V*.  w.  **,  r,  -  «.•  <-  V  V  v.v.v.%  .N  . 


Setting  equations  (13)  and  (14)  equal  to  each  other  and  using  equation  (8)  results 
in  an  expression  for  the  equivalent  mass  of  the  bimorph. 


M  -  0.242  Ma 


where 


Ma  «  pLWT  .  (15) 

With  the  density  of  the  bimorph  given  in  Table  2,  the  effective  mass  of  the 

-4 

bimorph  is  calculated  to  be  M  ■  3.8  x  10  kg. 

Evaluating  the  rystam  damping  is  not  as  easy.  Experience  has  shown  that  it 
is  fairly  small  relative  to  the  other  terms.  Experimental  methods  will  be  used 
to  determine  the  damping  coefficient  and  the  effective  bimorph  mass,  M.  If 
equation  (5)  is  solved  for  »  0.0  and  FX(S)  *  0.0,  the  resulting  characteristic 
polynomial  Is  given  by 

M  s2  +  B  $  +  l  *  0  .  (17) 


This  is  of  the  form 


where 


-Ls2  +  ^-S  +  l=  0  , 
WN2  WN 


WN  /  W 


the  natural  frequency  and  the  damping  ratio  can  be  obtained  from  the  experimentally 
measured  step  response  of  the  system.  Since  K$  is  known,  then  both  the  mass,  M, 


arid  the  damping  coefficient,  B,  can  be  determined. 

Figure  9  is  the  step  response  of  the  system  with  RA  =  0.0  and  Fx(t)  *  0.0. 

From  this  figure  the  damped  natural  frequency,  „>d,  is  seen  to  be  u>d  =  880  sec”1. 
Also,  using  logarithm  decrement  the  damping  ratio,  c,  can  be  shown  to  be  r,  =  0.005. 
The  natural  frequency  can  be  determined  from  the  following  relationship. 


V' '  -  <■' 


Figure  9.  Measured  System  Step  Response 

Ra  =  0.  OM.; 

F  (t)  =  0. ON 
x  ' 


\ 
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Since  the  damping  ratio  is  so  small,  the  natural  frequency  is  essentially  the 
same  as  the  damped  natural  frequency  and  s  880  s"\  The  measured  effective 
mass,  M,  of  the  system  is  then  calculated  to  be  M  =  3.8  x  10’^  kg.  The  analytically 
calculated  mass  using  equation  (15)  is  M  *  3.8  x  10"4  kg,  which  agrees  with  the 
measured  value.  The  system  damping  is  also  calculated  to  be  B  =  3.3  x  10  N-sec/m. 

Equation  (2)  shows  that  the  piezoelectric  force  constant,  C^,  can  be  calculated 
by  dividing  the  product  of  the  spring  constant,  K$,  and  the  steady-state  dis¬ 
placement,  X,  by  the  applied  voltage,  V.  Also  the  steady-state  displacement,  for 

an  applied  voltage,  at  the  end  of  a  bimorph  with  a  series  electrical  connection, 

12 

is  given  by  . 


X 


"2d31V72 


(22) 


The  piezoelectric  force  constant,  Cj,  is  expressed  by 


(23) 


C,  .  2Ksd3]  £  .  (24) 

-4 

For  the  bender  of  Table  2,  is  then  calculated  to  be  =  -8.6  x  10  N/v. 

The  piezoelectric  force  constant  can  be  verified  experimentally.  Figure  9 

shows  that  for  an  applied  voltage  of  14  V,  the  steady-state  displacement  is  measured 

to  be  4.0  x  10"5  m.  Applying  equation  (23)  results  in  a  measured  piezoelectric 

-4 

force  constant  of  =  -8.3  x  10  N/V.  This  compares  favorably  with  the  predicted 
value  of  C^. 

The  piezoelectric  back  current  constant  can  be  calculated  from  a  knowledge 

of  the  piezoelectric  bimorph's  properties  as  a  generator.  The  amount  of  charge 

1 2 

produced  is  proportional  to  the  applied  force.  This  relationship  is  given  by 


q 


3  fl  h 
2  ^2  a3! 


(25) 


12 


Piezoceramic  Design  Note,  Gulton  Industries,  Metuchen,  New  Jersey,  (1978). 
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Equation  (1)  shows  that,  in  the  steady  state  with  %  -  X  =  0  and  the  bimorph  open 
circuited,  the  applied  force  is  given  by 

FX  =  K$X  +  C,V  .  (26) 

The  voltage  potential  across  the  plates  of  the  bimorph  can  be  calculated  using 
Equation  (25),  (28),  (29),  (30)  and  is  given  by 

v  =  \  St  g3i  ’  ^ 

where  the  various  parameters  are  defined  as 

q  =  CTV  ,  (28) 

CT  =  e  ^  f  (29) 

G  a  d31/s31  ‘  (30) 

Substituting  equations  (27)  and  (26)  into  equation  (25)  resuts  in  a  relationship 
for  the  charge  produced  for  a  given  displacement  of  the  end  of  the  bimorph. 


q  3  s  L 

X  ~  2  ,  3r  L  _  T 
1  "  2  C1  WT  g31  1 


L2 

.2  d31 


Since  L,  W,  and  T  change  very  little  with  the  applied  force,  the  piezoelectric  back 
current  constant  is  given  by 


£  =  1 
x  2 


3  r  L_  _  T 
"  1  C1  WT  931  T 


Z?  d31 


(32) 


Substituting  in  values  for  t.he  test  bimorph  parameters  results  in  a  piezoelectric 

•4  A-s 

back  current  constant  of  K  =  -7.1  x  10  . 

q  m 

The  electrical  parameters  can  also  be  determined  analytically.  The 
capacitance  of  the  bimorph,  Cy,  can  be  calculated  by  combing  equations  (29)  and  (30), 
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d31LW 

=  <33> 

The  capacitance  of  the  test  bimorph  is  calculated  to  be  Cy  =  2.7  x  10"8F.  For  the 
purpose  of  the  model,  the  resistance  of  the  bimorph,  Ry,  is  assumed  to  be  infinite. 
Ry  was  measured  to  be  greater  than  10^  ohms,  which  is  very  large  compared  to 
RA,  so  this  is  a  good  assumption. 

The  electrical  parameters,  Cy,  Ry,  and  R^,  can  be  verified  experimentally. 

The  values  for  the  supply  resistor,  R^,  and  the  ferroelectric  materials  leakage  re¬ 
sistance,  Ry,  can  be  measured  with  an  ohm  meter.  The  bimorph  capacitance,  Cy, 
can  be  measured  as  any  ordinary  capacitor  with  the  motion  of  the  bimorph  constrained. 

O 

The  measured  value  for  the  bimorph's  capacitance  is  found  to  be,  Cx  =  3.0  x  10  F. 

The  leakage  resistance,  Ry,  is  so  large  that  it  is  assumed  to  be  open  circuited. 

The  supply  resistor,  R^,  is  changed  to  vary  the  frequency  response  characteristics 
of  the  bimorph.  All  of  the  measured  and  the  analytically  predicted  system  parameters 
are  sunwarized  in  Table  3. 

TABLE  3.  MEASURED  AND  ANALYTICALLY  PREDICTED  ELECTRO-MECHANICAL  SYSTEM  PARAMETERS 


Parameter 


System  spring  constant,  K$ 

Effective  mass,  M 

System  damping,  B 

Piezoelectric  force 
constant,  C^ 

Piezoelectric  back 
current  constant,  Kq 

Bimorph  capacitance,  CT 


Measured  Value 


290  N/m - 

3.8xl0“4  kg 
3.3xl0“3  N/m 


-8.3x10“*  N/V 


3.0xl0'8  F 


Analytically  Predicted 
Value 


310  N/m 

3.8x  10"4  kg 

5xl0’3  N/m  (estimated) 


-8.6x10“*  N/V 


7.1xl0"4  A-s/m 
2.6xl0’8  F 


Equation  (5)  is  then  simulated  on  the  computer  for  Fx (S )  *  0.0  and  VA(S) 

^A 

*  -$  »  a  step  voltage  input.  The  step  response  of  the  system  is  examined 
analytically  predicted  system  parameters.  The  anlaytically  determined  computer 
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output  is  then  compared  with  experimentally  determined  data.  Comparisons  between 
the  experimental  data  and  the  analytical  data,  for  various  values  of  the  supply 
resistor,  are  shown  in  Figures  10,  11,  12,  and  13. 

The  close  match  between  the  analytical  data  and  the  experimental  data  suggests 
that  the  model  is  sufficiently  accurate  for  design  work  in  the  frequency  range 
less  than  the  first  mode  of  resonant  vibration.  This  model  will  be  applied  in 
Section  5  in  the  design  of  an  electropneumatic  converter. 

4.  DEVELOPMENT  OF  THE  FLUID  MODEL 

A  fluid  model  describing  the  output  pressure  difference  of  this  LPA  for  a 

•  | 

given  input  motion  of  the  end  of  the  piezoelectric  bimorph  is  needed  to  match 

i 

a  gain  block  with  a  flapper-nozzle  system.  This  is  accomplished  by  combining  an 

i 

equivalent  steady-state  resistive  model  for  the  nozzles  with  a  modfel  for  the  LPA. 

•  i 

A  set  of  static  design  equations  for  the  LPA  is  then  used  to  design  the  LPA  gain 

i 

block.  Finally  the  gain  block  is  checked  for  speed  of  response.  In  this  manner 

i 

an  LPA  system  can  be  matched  to  the  nozzles.  j 

i 

l 

The  effect  of  static  and  dynamic  fluid  flow  forces  on  the  dynamic  response  af 

I 

the  bimorph  is  also  important.  Increasing  the  supply  pressure  to  the  bimorph 
chamber  causes  flow  forces  to  be  exerted  on  the  bimorph.  Therefore,  the  bimorph 
supply  pressure  should  be  reduced  to  a  level  that  maintains  system  stability. 

This  problem  ib  approached  analytically  with  a  control  volume  analysis  of  the 
flapper-nozzle  region. 

The  goals  of  the  fluid  model  are  twofold.  The  first  goal  is  to  demonstrate 
that  higher  supply  pressures  to  the  bimorph  will  result  in  a  potentially  unstable 
dynamic  response.  This  will  occur  because  of  the  static  and  dynamic  fluid  flow 
forces,  which  increase  with  higher  bimorph  supply  pressures.  The  other  goal 
will  be  to  design  an  amplifier  system  capable  of  a  given  output  pressure  and  flow 
with  a  given  system  banuwidth. 


infers  x  11) 


meters 


(b)  Experimental 

1  Large  division  horizontal  =  1.0 
1  Large  division  vertical  =  9.1  x 


Figure  12.  System  Step  Responses:  (a)  Simulated  and  (b)  E 
R.  *  0.260M.;,  F  ft)  =  0.0 
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The  flow  forces  exerted  on  the  bimorph  can  be  calculated  using  the  momentum 
equation  applied  to  the  control  volume  comprising  the  air  between  thr  bimorph  and 
the  nozzle.  Figure  14  illustrates  this  control  volume  along  with  the  necessary 
dimensions.  Applying  the  momentum  equation  for  the  forces  on  the  control  volume, 
including  the  force,  F,  exerted  on  the  control  volume  by  the  bender,  results  in 


-F  -  (P,  -  P,)  nr2  =  --r  i  oV  dV-  +  |  ,V  (V  •  n)  dA  .  (34) 

c  1  J c . V .  j c . s . 

This  vector  equation  is  then  rewritten  for  the  force,  F^,  exerted  by  the 
fluid  on  the  bender  in  the  direction  of  the  motion  of  the  end  of  the  bimorph. 

F  ■  Fx  ’  <P1  -  p2>^2  *  ft(*  UlLbl  -  ¥b?l  *  f+  ¥l  '  ¥2]  '  (35) 

Neglecting  compressibility,  neglecting  the  rate  of  change  of  volume  of  the  control 
volume  due  to  the  motion  of  the  end  of  the  bimorph,  and  applying  the  principle  of 
continuity  to  the  control  volum!e  of  Figure  14  yields 


1 

Wli  =  ko  =  W1  * 


W2i  =  ^20  “  W2  ' 

Imcompressible  flow  will  bq  assumed  because  of  the  low  operating  supply 

pressures  and  low  velocities.  Using  the  incompressible  flow  equations  from 
13 

Anderson  results  in 


=  Ce  . -r(u  +  X)  /2(PSi-P1)/r  , 

(38) 

=  Ce  .  -r(u  -  X)  /2(PS1  -  P2)/r  , 

(39) 

=  /2CPj~~Pyjr  ,  and 

(40) 

■  'Wv¥;  • 

(41) 

Substituting  in  for  each  component  in  the  momentum  equation  results  in 


JB.W.  Anderson,  The  Analysis  and  Design  of  Pneumatic  Systems,  John  Wiley  and 
Sons,  Inc.,  New  York,  (1976). 


I 


•-  V>  V*  V  T 


vr?-v.v^.vr^.'  •*"*  ^  *v  ^  t 
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W1V]  =  2Ce  -r(u  +  X)(P$1  -  P})  , 
W2V2  “  2Ce  ’r(u  ”  X)(PS1  "  P2}  * 


wiLbi s  ce  oLbrr(u  +  x)  /2{Psi  ~rjW:  ’  and 

W2Lb2  *  Ce£Lb2'r<u  -  X)  /2<PS1  -  P2)/;  • 


where 


Lbi  *  4,  *  “  * x  • 


•b2 


Lb  +  u  •  X 


(42) 

(43) 

(44) 

(45) 

(46) 

(47) 


Taking  the  partial  derivative  with  respect  to  time  of  equations  (44)  and  (45) 
results  in 

It  ‘ulLbl)  *  (ce  £Lbl  =r  ''2<PS1  -  V/;)  *  * 


(48) 


+  (Ce  ctr  (u  +  X)  /2(PS1  -  P,)/.-)  *  + 


-  (Ce  cLbl  -r(u  +  X)  ( 


1/2  <PS1  -pi> 


It,  . 


ft  <“2Lb2>  ■  -  <Ce  £Lb2  -^2(PS,  -  I •!>/■-)  *  * 

-  (Ce  c-r(u  -  X)  /2(PS1  -  P i ) / >' )  *  + 


(49) 


(C  .L  -r(u  -  X)  (———))  t2 
e  b2  /2.(Ps“Pj  2 


since 


and 


L.  >>  u  +  X 
o 


XSX 


max 


Lh  »  u  *  XI 


(50) 

(51) 


'X=X  . 
min 


Therefore,  the  second  position  rate  term  in  equation  (48)  and  (49)  is  neglected. 


\ 


/: 


J. 
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McCloy  and  Martin  state  that  the  pressure  rate  term  is  small  in  comparison  with  the 

14 

remaining  term,  so  it  is  also  neglected  . 

Making  the  final  substitution  into  the  momentum  equation  results  in 

FX  =  (P,  -  P2)  +.  Ce  Lb  ^r(/2(Psl  -  FjT/p  +  2(PS1  -  P2)/p]  *  + 

+  2  Ce  |r  +(u  +  X ) ( P $ -j  -  P,)  -  (u  -  X)(PS1  -  P2)  (52) 

The  pressures  P^  and  P^  must  be  known  before  the  momentum  equation  can  be 
solved.  A  steady-state  fluid  model  is  first  analyzed  to  determine  these  two 
pressures.  Figure  15  is  an  equivalent  resistive  circuit  for  the  nozzles  and  the 
LPA  input. 

Using  the  principle  of  continuity  at  the  P^  node  and  the  P^  node,  respectively, 
results  in 

“l  ■  “e  *  UA1  •  <53> 

w2  ■  %  •  V  •  (54) 

2 

Using  the  equations  developed  by  Manion  and  Drzewiecki  for  the  LPA  input 
results  in  an  expression  for  the  jet  entrainment  flow, 

We  =  iQe  ,  (55) - 

where 

1  CnBr)l/3  i 

Qe  -  0.021  CdC)}  +  -  1  Qs  .  (56) 

and  ~ 

QS  *  Cdhbs  /2P7/7  ,  (57) 

Nr  *  (bs/a)  72(Ps  -  Pv)/p  .  (58) 

14 . . .  “  —  '  ’ 

0.  McCloy  and  H.R.  Martin,  The  Control  of  Fluid  Power,  John  Wiley  and  Sons,  Inc. 

New  York  (1973). 

F.M.  Manion  and  T.M.  Drzewiecki,  Analytical  Design  of  Laminar  Proportional 
Amplifiers,  Proc.  HDL  State-of-the-Art  Symposium,  Vol .  1,  Adelphi,  MO  (Oct.  1974). 
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Solving  for  the  LPA  input  channel  resistance^  results  in 


24  X  C  . 
c  a  a 


~  1  +  ~  +  0.35  R.  , 


where 


c  V  nr  Bc 


r s  =  r  •  (60) 

ws 

2 

The  LPA  vent  resistance  ,  R v,  has  been  found  to  equal 

Ry  =  3.0  R$  .  (61) 

Returning  to  Figure  15  and  solving  for  the  flow  escaping  from  the  vents,  on 
side  1,  results  in 

rpn  ‘  paI 

“Al  ‘  — "j  <!  (62) 

where 

•Vi  s  pi  -  r  “i  •  <63) 

Substituting  equation  (63)  into  equation  (62)  and  solving  for  the  vent  flow 
using  gauge  pressures  results  in 


Similarly, 


Al  ~  R 


A  2  R 


f(Y 

R 

fU 

P  - 

No 

R 

Vj 

vj 

( p 

2 

"* 

fy 

y 

Rv 

Sutbsti tuting  equation  (64)  and  equation  (66)  into  equation  (53)  and  equation 
(54),  respectively,  results  in 

f  Rc]  [Pi] 

W1  !  1  +  -  r1|  -  w  =  0  .  (661 


“2  n  ‘  «'  - 


V 

i  ■ 

R 

i 

v 

^ . 

i 

,pi 

“ 

R 

v; 

■  ■  -  «e  *  0  . 


I.  -  W  =0 
e 


Equations  (52),  (66),  and  (67)  describe  the  fluid  flow  forces,  F^,  exerted  on 
the  bender  by  the  nozzles  and  the  pressures,  P^  and  P^,  downstream  of  the  nozzles. 
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These  equations  are  a  set  of  nonlinear  differential  equations  which  can  be  linearized 
about  an  operating  point  and  then  applied  to  the  mechanical  model.  Rewriting 
equations  (52),  (66),  and  (67)  by  substituting  for  the  control  flows  and  separating 
the  inputs,  Xi  X,  from  the  outputs,  P^,  P2,  Fx,  results  in 


Wj  (Pr  X)  »  A,  (u  +  X)  -  Pi  -  A,P1  -  A,  -  0  , 


SI  rl  21 


(68) 


»'2  (P2,  X)  •  A]  (u  -  X)  /PS1  -  P2  -  A2P2  -  A3  *  0  ,  (69) 

Fx  (Pr  P2,  X,  *)  3  A4  (P}  -  P2)  +  Asi  (/Psl  -  Pj  +  /PS1-  P2)  + 


+  Ag(+  (u  +  X)(P$l  -  P,)  -  (u  -  X)(P$1  -  P2)), 


where 


A]  3  Ce  pnr 


A2  *  p/Rv 


A3  =  We 
A4  3  nr* 


f,  *  y 


S2/p 


A5  =  Ce  pLb  *r  ^ 

A,  =  2C„7ir 
— o - e 


(70) 

(71) 

(72) 

(73) 

(74) 

(75) 

(76) 


Linearizing  equations  (68),  (69),  and  (70)  about  >j  *  P2  3  0,  X  *  0,  X  =  0 
defines  an  operating  point  where  the  bimorph  is  stationary  and  undeflected  and  the 
laminar  jet  is  centered.  Solving  for  this  point  results  in 


Alu  l^s7"  P1o  "  A2P1o  *  A3  =  0 


Alu  '^Sl  '  P2o  "  A2P2o  -  A3  =  0 

Fx  =  0  • 


(77) 

(78) 

(79) 


Because  of  symmetry  at  the  operating  point  which  can  be  noted  in  equations  (77) 
and  (78),  the  operating  point  bias  pressures,  P1(J  and  P2q,  are  equal. 
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(89) 


\ 

\ 
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Similarly,  evaluating  equation  (83)  and  solving  for  the  pressure  at  nozzle  2 
results  in 


P2  = 


2A,uP2  -  2A,P  +  A,uP„  -  2A,P2X 
l _ 3 _ 1  0 _ 1 

A-|U  +  2A2P 


where 


pa/psTVTo 


(90) 


(91) 


Evaluating  equation  (84)  results  in  a  linearized  expression  for  the  system 
fluid  flow-forces. 

FX  *  (A4  -  Agu)(P1  -  P2)  +  2Ag  P2  X  +  2A5  P  %  .  (92) 

Substituting  equations  (89)  and  (90)  into  equation  (92)  results  in  a  final 
expression  for  the  linearized  system  fluid  flow  forces  as  a  function  of  the  position 


and  the  velocity  of  the  end  of  the  bimorph. 


FX  * 


4A1  Pd  (A4  -  A6u) 


A-jU  +  2A^P 


+  2AgP* 


►  * 

X  + 

- 

2a5P 

Jt  . 


(93) 


Taking  the  Laplace  transform  of  equation  (93),  setting  the  initial  conditions 
equal  to  zero,  and  solving  for  FX(S)  gives 

FX(S)  *  X(S)(C0  +  C2S)  .  (94) 

The  constant  Cg  is  the  fluid  system  spring  constant  and  C2  is  the  fluid  system 
damping  coefficient. 


Co  * 


4AjP 


2  f 


A'<  "  + 


A^u  +  2A2P 


(95) 


^2  *  ^5P 


(96) 


Substituting  equation  (94)  into  equation  (5)  and  solving  for  X(S)/V^(S) 


gives 


X(S)  . 


-v 

1 

o 

o 

o 

4* 

1  c 

L6 

s3  ♦ 

C7  ‘  C2C3 

2 

S^  + 

C8  -  C0C3  -  C2C4 

S  +  1 

1  "  C0C4 

1  ‘  C0C4 

,  (97) 


where 
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ctr 

c  =  1 

C3  K$  ‘ 

(98) 

* 

it 

(99) 

rt  c, 

cs  ks1ra  +  rtJ  • 

(100) 

ctmr 

r  =  1 

6  Ks  ’ 

(101) 

CtBR  +  M 

C7  Ks 

(102) 

„  <CTKs  *  ClVR  +  B 

L8  K$ 

( 1 03} 

It  has  been  experimentally  determined  that  increasing  the  supply  pressure  to 
the  bimorph  has  a  destabilizing  effect  on  the  system's  dynamic  response.  Equation 
(97)  will  be  used  in  Section  5  with  specific  LPA  and  bimorph  parameters  to  demon¬ 
strate  that  this  is  a  result  o*  the  static  and  dynamic  fluid  flow- forces  exerted 
on  the  bimorph  by  the  nozzle. 

Equations  (89)  and  (90)  can  be  used  to  predict  the  pressures,  and  P £,  which 
are  output  from  the  flapper-nozzle  and  are  the  input  co  the  LPA.  Once  P-j  a^ 

?2  are  known  as  a  function  of  the  flapper-nozzle  configuration,  a  set  of  LPA  design 
equations  can  be  utilized  to  amplify  this  pressure  signal  to  a  usable  level.  The 
LPA  design  equations  are  a  set  of  steady-state  design  equations.  The  LPA  design 
method  used  to  devise  an  LPA  gain  block  is  to  design  it  statically  and  then  go  back 
and  check  if  it  is  fast  enough  for  the  intended  purpose. 
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The  steady  state  LPA  design  equations  for  pneumatic  application  are  given  by 


Reynolds 

number, 

dimensionless 

Nr  =  (bs/v)/2(p$  -  Pv  /p 

(104) 

Modified 

Reynolds 

number 

dimensionless 

_  nr 

Nr  ~  iT?  ’ 

2  1  +  o  J 

(105) 

Supply 
pressure, 
mm  Hg 

p  a  »R  ’2V2V5  2 

Ks  "  H  120  0  Bs 

(106) 

Supply  flow, 

LPM 

rNA  1  b  ' 

Qs  '  0  3  [l20j  (o. 5  * 

(107) 

Supply  nozzle  fluid 
resistance, 
mm  Hg/LPM 

fw* 

NR  1  2  0.5  2 

KS  U  [l  20j  [oj  bs 

(108) 

Amplifier  input 
resistance, 
mm  Hg/LPM 

Ri  =  0.75  Rs 

(109) 

Amplifier  output 
resistance, 
mm  Hg/LPM 

Rq  =  0.50  R$ 

(no) 

Aspect  radio 
dimensionless 

o  -  h5/bs 

(HD 

Gain, 

dimensionless 

gp  ■  gpo  lyf1  * 

(112) 

These  equations  will  produce  results  in  the  units  shown  if  bs  is  given  in  mm. 

To  maximize  the  gain  of  the  amplifiers,  the  modified  Reynolds  number,  NR,  is  chosen 
as  large  as  possible  while  keeping  the  LPA  in  a  laminar  region.  The  maximum  modified 
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Reynods  number,  N^,  for  laminar  flow  in  the  LPA  has  been  determined  by  Manion 
to  be  120.  For  maximum  gain,  each  stage  of  amplification  will  then  operate  at  a 
supply  pressure  that  will  result  in  a  modified  Reynolds  number  of  120. 

The  dynamic  range  of  an  amplifier  is  defined  as  the  ratio  of  the  largest 
input  to  the  smallest  input  to  the  amplifier  that  gives  uniform  gain.  For 
maximum  dynamic  range,  each  stage  of  the  LPA  gain  block  should  saturate  at  the 
same  time.  This  means  that  for  a  given  pressure  difference  input,  each  jet  in 
each  stage  should  sweep  the  same  angle.  This  will  prevent  premature  saturation 
of  the  entire  gain  block  due  to  saturation  of  an  intermediate  stage.  This  can  be 
accomplished,  as  shown  below,  utlizing  the  LPA  design  equations^. 

If  the  gain  block  is  considered  to  be  comprised  of  N  stages  of  amplification, 
then  the  output  pressure  of  the  N  stage  equals  the  input  pressure  to  that  stage 
multiplied  by  its  gain.  Therefore, 


pressure 


P0  (N  -  1)  -  P0(N)/Gp  . 

Moreover,  since  the  supply  pressure  is  linearly  related  to  the  output 
16 


013) 


Ps  (N  -  1)  =  P$  (N)/Gp 


Equation  (106)  gives, 

i  NL  (N-l )  2  •.  c  q  c  2 

Ps(N  \  =  4  (~T20  )  (o(N-lp  (b  (N-l)* 


and 


PS(N)  =  4 


120 


(  1 
^W) 


1; 


bfwi  *  Vs  <*-»• 


015) 


016) 


If  all  N  stages  of  amplification  are  assumed  to  operate  at  the  same  modified 
Reynolds  number  and  each  stage  has  the  same  nozzle  width,  then  the  optimum  aspect 
ratio  for  maximum  dynamic  range  is  derived  by  combining  equations  (115)  and  (116). 
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This  is  given  by 

f$rr°-3  •  <m> 

where 

6p-10. 

These  are  generally  good  assumptions  because  the  modified  Reynolds  number,  N^,  is 
usually  fixed  as  high  as  possible,  while  still  remaining  in  the  laminar  region, 
to  maximize  the  amplifier’s  gain.  Also  the  nozzle  width  would  remain  constant  if 
identical  amplifier  laminations  were  used  for  each  stage.  Therefore,  the  easiest 
design  parameter  to  vary  to  satisfy  equations  (114),  (115)  and  (116)  is  the 
aspect  ratio  since  this  can  be  adjusted  by  adding  or  taking  away  amplifier 
laminations. 

When  applying  equation  (117)  It  Is  important  to  note  that  there  are  limitations 
on  the  size  of  the  aspect  ratio.  Orzewlecki16  has  determined  that  amplifiers  with 
an  aspect  ratio  of  less  than  0.3  have  reduced  gain  due  to  viscous  losses.  There¬ 
fore,  allowing  o(N)  *  0.3  as  the  smallest  allowable  aspect  ratio  for  the  final 
stage  of  amplification,  the  optimum  aspect  ratios  of  a  gain  block  can  then  be 
calculated.  This  results  in  aspect  ratios  of  o(N)  *  0.3,  o(N-l)  *  1.0,  o(N-2)  = 

3.3,  o(N-3)  *  11.1,  etc.  Aspect  ratios  for  gain  blocks  of  more  than  three  stages 
are  too  large  for  practical  purposes.  Therefore,  a  gain  block  of  three  stages  is 
the  largest  practical  size  which  can  be  built  if  the  same  amplifier  laminations 
are  to  be  used  for  each  stage  and  the  system  is  to  have  maximum  dynamic  range. 

If  more  gain  is  required  of  a  system  than  can  be  achieved  from  just  three 
stages  of  amplification,  the  supply  nozzle  width  can  also  be  varied  to  design  a 
gain  block  with  maximum  dynamic  range.  Each  supply  nozzle  width  would  have  three 
stages  of  amplification  associated  with  it.  These  three  internal  stages  would 
have  aspect  ratios  calculated  using  equation  (117).  If  an  aspect  ratio  becomes 


6T.M.  Drzewiecki ,  A  Fluidic  Audio  Intercan,  20th  Anniversary  of  Fluidics 
Symposium,  The  ftnerican  Society  of  Mechanical  Engineers,  (1981). 
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too  large  or  too  small  for  practical  purposes,  equation  (116)  could  be  used  to 
calculate  a  new  optimum  supply  nozzle  width,  and  the  staging  could  continue. 

An  example  of  this  would  be  if  four  or  more  stages  of  amplification:  were 
necessary  for  a  given  application.  This  example  would  require  the  designer  to 
pick  an  initial  supply  nozzle  width,  bs,.and  use*  that-value  for  the  first  three 
stages  of  amplification  with  0(N)  =  .0.3,.  o(N-l)  =  1.0,  o(N-2)  «  3.3'"  To  maximize 
the  dynamic  range  of  the  next  set  of  three  amplifier  stages,  the  aspect  ratio  of 
the  (N-3)  stage  is  chosen  to  be  a(N-3)  **  0.3.  -Applying  equation  (116)  to  this 
results  in  an  expression  for  the  supply  nozzle  width  ratio  for  this  next  stage. 


-•-'"V 


b$(m) 

bsfm-T) 


=  0.03 


(118) 


The  number  of  three-stage  amplifier  sets  under  consideration  in  equation  (118)  is 
designed  by  m.  Stages  (N-4)  and  (N-5)  can  then  be  designed  with  a(N-4)  ,»  1.0,  an 
o(N-5)  *  3.3.  This  process  can  continue  in  this  manner  until  enough  amplification 
has  been  realized.  The  disadvantage  with  this  approach  is  that  the  physical 
size  of  the  amplifiers  becomes  large  quickly,  making  a  large  gain-block  prohibitive. 
Therefore,  other  techniques  for  designing  a  gain  block  with  more  than  three  stages 
is  necessary. 

Large  gain  blocks  can  be  fashioned  in  a  manner  which  places  less  of  a 

* 

penalty  dynamic  range  and  stresses  other  factors  which  may  be  important  to  the 
system  design.  The  two  best  methods  for  this  type  of  design  are  a  self-staging 
technique  and  a  common-sense  approach  which  uses  a  knowledge  of  the  input  signal 
size  and  avoids  a  gain  block  that  saturates  for  those  inputs. 

The  self-staged  amplifier  gain  block  is  a  set  of  identical  amplifier  stages 
cascaded  together.  The  advantages  of  self-staging  are  simplicity  of  design  and 
the  ability  for  increased  gain.  There  is  only  one  pneumatic  supply  in  a  self- 
staged  gain  block  and,  therefore,  no  need  for  dropping  resistors  to  supply  each 


separate  stage.  The  resolution  of  the  LPA  is  sufficiently  high  so  that  an 
amplifier  may  be  self-staged  many  times  before  any  significant  reduction  is  output 
resoltuion  is  detectable.  The  disadvantages  of  a  self-staged  amplifier  are  a 
lower  dynamic  range,  since  the  first  jet  in  the  gain  block  will  sweep  a  much 
smaller  angle  than  the  last  jet,  a  smaller  usable  frequency  range,  and  lower  gain 
per  stage  due  to  a  poorer  input-output  impedance  match. 

The  final  approach  to  LPA  gain  block  design  is  for  special  purpose  amplification. 
Occasionally,  achieving  maximum  dynamic  range  is  not  as  important  as  another  de¬ 
sign  criterion.  The  designer  may  desire  to  operate  each  stage  at  as  high  a  supply 
velocity  as  possible  to  minimize  transport  delay.  The  designer  may  also  choose  to 
design  LPA  stages  with  an  improved  input-output  impedance  match.  This  may  be 
necessary  to  maintain  reasonable  gain  within  a  gain  block  with  many  parallel 
stages.  In  any  case,  the  approach  used  to  devise  a  special  purpose  amplifier  is  to 
identify  the  characteristics  which  are  most  important  in  the  design.  The  gain 
block  is  then  designed  analytically  and  the  trade-offs  for  the  special  purpose 
application  are  examined.  If  the  trade-offs  are  acceptable,  the  design  is  used. 

The  LPA  gain  block  system  is  designed  using  one  of  the  approaches  described 
above.  Once  an  approach  has  been  selected,  the  next  points  of  interest  are 
static  gain  and  frequency  response.  The  static  gain  of  one  stage  of  amplification 
is  given  by  equation  (112).  The  total  gain  of  the  gain  block  is  the  product  of 
the  gains  of  each  stage.  The  output  pressure  difference  of  the  gain  block  is  then 
the  input  pressure  difference  times  the  gain. 

The  input  pressure  difference,  for  a  given  bimorph  displacement  to  the  gain 
block,  can  be  derived  by  combining  equations  (89)  and  (90).  This  results  in  an 
equation  describing  the  input  pressure  difference  to  the  gain  block  for  a  steady- 
state  position  of  the  end  of  the  bimorph. 


rVV  ' 


1  ■  2 
X 


(119) 


The  pressure  difference  output  of  the  nozzle  for  an  input  deflection  of  the 
end  of  the  bimorph,  expressed  by  equation  (119),  can  be  used  to  design  the  first 
stage  of  amplification.  The  linear  pressure  recovery  fo  an  LPA  stage  is  given 
by  Drzwiecki  to  be  ^ 


P01  ‘  P02 


=  0.70  P. 


where 


P01  "  P02  *  GP 


pii  -  pi2] 

) 


(120) 


021) 


It  is  important  that  the  maximum  input  pressure  difference  is  matched  to  the 
amplifier  gain  block  to  avoid  saturation  in  any  of  the  stages.  The  supply  pressure, 
P$,  to  an  amplifier  stage  is  given  by  equation  (106).  Equation  (106)  also  shows 
that  the  supply  pressure,  Ps,  is  inversely  proportional  to  the  square  of  the 
aspect  ratio,  the  most  easily  varied  design  parameter.  Therefore,  for  a  given 
input  pressure  difference,  the  first-stage  aspect  ratio  may  be  varied  to  achieve 
a  matched  first-stage  amplifier,  where  a  maximum  sweep  of  the  first  stage  jet  is 
realized  and  the  amplifier  gain  block  does  not  saturate - 


\  \  ... 


Each  additional  stage  is  designed  using  an  optimum  dynamic  range  staging 
technique,  self-staging  technique,  or  a  special  purpose  technique.  The  gain  of 
the  system  is  then  calculated  using  equation  (112)  for  each  stage. 

Equation  (112)  states  that  the  single-stage  gain  of  an  amplifier  is  a  function 
of  the  blocked  load  gain,  Gpg,  the  amplifier's  output  impedance,  Rq,  and  the  impedance 
of  the  load,  .  In  a  cascaded  network  the  load  impedance  of  the  (N-l)  stage  is 
equal  to  the  input  impedance  of  the  (N)  stage.  The  blocked  load  gain,  Gpg,  of 
standard  LPA  amplifiers  varies  between  9  and  10  per  stage.  With  this  information. 


T.M.  Drzewiecki,  A  Fluidic  Voice  Communication  System  and  Data  Link,  Ph.D. 
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one  may  use  equations  (108),  (109),  (110),  and  (112)  to  calculate  the  gain  of  an 
LPA  gain  block. 

The  final  consideration  of  the  amplifier  system  is  its  frequency  response. 

The  phase  response  component  of  the  system  frequency  response  is  generally  dominated 

2  16  17  18 

by  signal  transport  lag  between  stages  The  phase  lag  is  derived  by 

calculating  the  time  for  propagation  of  the  signal  through  the  amplifier  system. 

2  16  17  18 

The  magnitude  response  of  the  system  has  been  determined  by  Drzewiecki  * 
and  a  guide  for  the  upper  limit  of  the  usable  frequency  range  for  uniform  gain  of 
Z-identical,  self-staged  amplifiers  is  given  by 


f  ,  o.03  ^  Z‘2/3 


(122) 


Equation  (122)  demonstrates  that  self-staging  an  amplifier,  with  Z  >  1,  will  result 

In  a  lower  usable  frequency  range  than  the  single-stage  amplifier  would.  The 

technique  for  calculating  the  usable  frequency  range  for  an  LPA  gain  block  is  to 

calculate  the  usable  frequency  range  for  each  stage  of  amplification.  The  dominant 

stage  would  be  the  stage  with  the  lowest  usable  frequency  range;  this  frequency 

value  can  be  used  as  a  guide  to  the  usable  range  of  the  entire  gain  block.  A 

more  exact  representation  of  the  dynamic  behavior  of  the  LPA  can  be  made  by  examining 

2  161718 

the  work  of  Manion  and  Drzewiecki  and  Drzewiecki  *  *  and  applying  one  of  the 

models  presented  in  those  papers  to  each  stage  of  the  gain  block.  Once  the 
gain  as  a  function  of  frequency  for  each  LPA  stage  is  determined,  the  overall 
LPA  gain  block  frequency  response  can  be  calculated.  This  is  the  product  of  the 
gains  of  each  amplifier  stage  as  a  function  of  frequency. 
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The  magnitude  response  may  also  exhibit  a  pronounced  resonant  peak  at  a  fre¬ 
quency  less  than  the  frequency  calculated  by  equation  (122).  This  phenomenon 

17  i  o 

is  explained  by  Drzewiecki  ’  as  a  result  of  internal  acoustic  feedback.  The 
mechanism  of  this  feedback  is  described  as  follows.  As  the  laminar  jet  oscillates 
within  the  LPA  geometry,  there  is  a  spilling  of  flow  to  one  side  and  then  to  the 
other.  The  flow  spilled  to  the  vents  is  radiated  acoustically  to  the  control 
region  where  it  may  combine  with  the  control  signal  in  the  form  of  positive  feed¬ 
back.  When  the  control  input  frequency  equals  the  acoustic  feedback  frequency, 
one  would  expect  that  there  would  be  an  increase  in  the  amplifier's  gain.  To 
determine  the  feedback  frequency,  one  must  first  note  that  the  jet  deflection 
travels  at  one-half  the  average  jet  particle  velocity.  Therefore  the  signal  is 
delayed  by  twice  the  jet  transport  time,  T',  and  the  following  equation  can  be 
written  for  the  period,  :,  of  the  acoustically  caused  oscillation.  Note  that 
travel  time  to  the  splitter  is  T1,  but  the  feedback  signal  returns  to  the  control 
region  at  the  speed  of  sound.  One  complete  cycle  occurs  in  T  s,  where  a  delay, 

T',  is  experienced  for  a  deflection  in  one  direction  and  a  delay,  T',  is 
experienced  for  a  deflection  in  the  other  direction. 

\:  =  2V  ,  (123) 

where 

r  -  xsp/ufl  ,  024) 

Ufl  *  Cd.»?r  .  (125) 

The  frequency  of  oscillation  due  to  internal  feedback  can  be  calculated  from  one 
over  the  period  and  is  given  by 

^T.M.  Drzewiecki,  A  Fluidic  Voice  Communication  System  and  Oata  Link,  Ph.D.  Thesis 
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(126) 


Ff  =  1/T  . 

Equation  (126)  can  then  be  used  to  calculate  the  frequency  at  which  resonant  peaks 
from  the  acoustical  feedback  do  occur.  Equation  (126)  will  be  applied  in  Section  5 
to  calculate  acoustically  caused  resonant  peaks  in  the  frequency  response  of  the 
working  prototype. 

The  fluidic  modeling  in  this  section  has  two  main  purposes.  The  flow  force 
modeling  is  important  because  it  will  be  used  to  show  that  increasing  the  bimorph 
supply  pressure  has  a  destabilizing  effect  on  its  dynamic  response.  The  equivalent 
resistive  fluid  circuit  analysis  coupled  with  the  LPA  design  equations  allow  the 
design  of  an  LPA  gain  block  with  output  characteristics  that  match  the  design 
intentions.  Using  the  results  of  this  section  and  Section  3,  one  can  design  and 
build  an  electropneumatic  converter.  The  analysis,  construction,  and  testing  of 
a  prototype  design  will  be  discussed  in  the  next  section. 

5.  DESIGN  AND  TESTING  OF  A  VALVE  SYSTEM 

The  equations  developed  in  Sectio..„  3  and  4  can  now  be  used  to  design  a  valve 
system.  A  summary  of  the  important  physical  proper  ies  of  the  prototype  valve 
system  is  provided  in  Appendix  A  for  reference.  Once  a  valve  system  has  been 
designed  analytically,  a  working  prototype  is  built.  The  working  prototype  is 
then  tested  to  demonstrate  the  validity  of  the  design  equations.  The  experience 
gained  through  this  process  results  in  design  changes  that  will  give  a  better 
system  response.  These  design  changes  are  discussed  in  Section  6. 

5.1  Analytical  Design 

The  first  step  in  designing  a  valve  system  is  to  identify  what  is 
important  in  the  intended  design.  The  static  and  dynamic  response  of  the  system 
is  a  function  of  the  LPA  gain  block  and  the  flapper-nozzle  component.  The  re¬ 
quired  response  of  the  system  can  be  determined  from  a  knowledge  of  the  intended 
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system  application.  Once  the  valve's  input  and  output  requirements  are  identified, 
a  system  can  be  designed. 

The  most  important  asp  .  of  building  the  first  working  prototype  of  a  new 
design  is  to  demonstrate  that  the  idea  works.  It  is  also  desirable  to  show  that 
an  organized  method  exists  for  improving  that  first  design;  in  this  case  that 
means  verifying  the  design  equation  of  Sections  3  and  4.  Once  seme  ...  fidence  in 
the  design  equations  has  been  established,  these  equations  may  be  used  for  de¬ 
signing  new  prototypes  for  specific  apDl ications. 

Since  the  LPA  defines  most  of  the  output  characteristics  of  the  system,  de¬ 
signing  the  LPA  gain  block  is  usually  the  initial  analytical  step  in  designing  a 
valve  system.  Before  the  LPA  gain  block  can  be  designed  however,  certain 
flapper-nozzle  parameters  must  be  assumed.  This  is  because  there  are  also  certain 
design  limitations  placed  on  the  bimorph.  The  diameter  of  the  nozzles  should  be 
as  small  as  possible  to  minimize  the  fluid  flow  forces.  The  diameter  of  the 
nozzles  will  then  limit  the  minimum  width  of  the  bimorph.  The  generally  accepted 

industry  standard  for  the  smallest  diameter  nozzle,  to  avoid  clogging,  is  D  * 

-4 

8.6  x  10  m.  Using  highly  filtered  air  this  dimension  could  be  reduced  further. 

To  insure  a  good  overlap  of  the  bimorph  and  the  nozzle,  the  bimorph  width  is 

_3 

then  chosen  to  be  W  =  1.4  x  10  m.  Also  the  bimorph  is  commercially  available 

-4 

in  a  single  standard  thickness.  That  thickness  is  given  to  be  T  *  4.9  x  10  m. 

With  the  width  and  the  thickness  of  the  bimorph  specified,  the  only  re¬ 
maining  bimorph  design  parameter  is  its  length,  L.  Equations  (13)  and  (16)  show 
that  the  first  mode  natural  frequency  of  the  bimorph  is  proportional  to  one  over 
its  length  squared.  Since  speed  of  response  is  an  important  design  criterion, 
it  is  desirable  to  reduce  the  length  of  the  bimorph  as  much  as  possible.  A 
shorter  bimorph,  though,  will  result  in  a  reduced  output  displacement  for  a  given 
applied  voltage.  The  smaller  displacement  will  result  in  a  reduced  pressure 
difference  at  the  nozzles.  This  will  require  additional  pneumatic  amplification 
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to  raise  the  pneumatic  signal  to  the  level  obtained  with  the  longer  bimorph. 

The  additional  pneumatic  amplification  can  be  delivered  with  more  LPA  stages. 

The  larger  gain  block  will  have  a  slower  frequency  response.  This  reduction  in 

frequency  response  is  caused  by  increased  signal  path  length  between  elements  and 

a  lower  bandwidth  found  in  low  supply  pressure  amplifiers  which  are  necessary 

to  amplify  the  reduced  output  of  the  nozzles.  There  is  an  optimum  obtained  by 

trading  off  bimorph  size  for  number  of  amplifier  stages. 

The  optimum  length  for  a  bimorph  is  therefore  a  function  of  the  gain  block 

used.  The  gain  block  describes  the  output  characteristics  of  valve  system  and 

will  therefore  vary  with  the  intended  application.  For  the  purpose  of  testing 

the  design  equations  and  building  a  working  prototype  of  the  system,  it  is  not 

necessary  to  specify  the  gain  block  first.  Therefore  a  bimorph  with  a  practical 

length,  resulting  in  a  moderately  high  frequency  response  and  a  moderately  high 

_2 

output  displacement,  is  chosen.  That  length  is  chosen  to  be,  L  *  1.5  x  10  m. 

The  natural  frequency  of  a  bimorph  of  this  length  is,  *  1050  Hz,  and  its  maximum 
output  displacement  for  a  15  V  p-p  input  is  -5.0  x  10*^'  m.  Once  the  validity 
of  the  design  equations  is  shown,  a  designer  could  optimize  the  length  of  a  bimorph 
for  a  specific  LPA  output  requirement. 

_ Once  the  bimorph  is  chosen *_±ts._specific  and  physical  material  properties 

are  defined.  These  properties  are  given  in  Table  4.  The  equations  developed  in 
Section  3  are  then  applied  to  analytically  predict  the  electromechanical  system 
parameters.  The  results  of  these  calculations  are  shown  in  Table  5. 

Now  that  the  bimorph  has  been  chosen  and  its  properties  are  identified,  an 
LPA  gain  block  can  be  designed  for  pneumatic  amplification  of  the  pressure  signal 
generated  downstream  of  the  nozzles.  A  standard  Corning  made  LPA  gain  block  is 
chosen  to  verify  the  analysis.  The  resulting  output  of  this  gain  block  is  used  to 
verify  the  design  equations  and  show  their  usefulness  for  developing  an  improved 
design  for  future  application.  The  purposes  of  this  prototype  design  and  working 


58 


TABLE  4.  PHYSICAL  AND  MATERIAL  PROPERTIES  OF  LEAD  ZIRCONATE  TITANATE  CERAMIC 
BTMORPH  USED  FOR  WORKING  PROTOTYPE 


Property 

Measured  Value 

Length,  L 

1.5  x  10'2  m 

Width,  W 

1.4  x  10"3  ffl 

Thickness,  T 

4.9  x  10‘4  m 

Nickel  plate  thickness,  tj 

1.52  x  10"6  a 

Lead  zirconate  titanate  thickness,  t2 

1.93  x  10‘4  m 

Brass  shim  thickness,  tg 

1.03  x  10‘4  m 

Density,  p 

7.5  x  103  kg/m3 

Piezoelectric  charge 
coefficient,  dg^ 

-1.8  x  10"10  m/V 

Piezoelectric  voltage 
coefficient,  gg^ 

i 

-1.1  x  ID'2  V-m/N 

Nickel  elastic  modulus,  Ej 

1.26  x  1011  N/m2 

Lead  zirconate  titante 

1.09  x  1011  N/m2 

elastic  modulus,  E2 

Brass  elastic  modulus.  Eg 

9.0  x  1010  N/m2 

TABLE  5.  ANALYTICALLY  PREDICTED  SYSTEM  PARAMETERS  USED  FOR  WORKING  PROTOTYPE 

Parameter 

Analytically  Predicted  Value 

System  spring  constant,  K$ 

832  N/m 

Effective  mass,  M 

1.9  x  10"5  Kg 

System  damping,  B 

5  x  10~3  N-s/m  (estimated) 

Piezoelectric  force 
constant,  C^ 

2.8  x  10"4  N/V 

Piezoelectric  back 
current  constant,  K^ 

-2.3  x  10'4  A-s/m 

Bimorph  capacitance,  C-j. 

1.1  x  10‘9  F 

model  are  to  demonstrate  that  the  device  does  work  and  that  the  design  equations 
do  predict  its  response. 

Therefore  an  amplifier  system  with  a  three-stage  preamplifier  and  a  single- 
stage  output  amplifier  is  chosen  for  the  design.  A  three-stage  preamplifier  is 
chosen  because  it  would  have  sufficiently  high  bandwidth  and  sensitivity  to  detect 
and  amplify  the  low  pressure  signal  coming  from  the  nozzles.  The  three  stages  of 
the  preamplifier  are  self-staged.  Each  stage  is  a  standard  LPA  design  with  its 
necessary  geometric  properties  given  in  Table  6.  The  jet  deflection  angle  of  the 
amplifier  is  then  designed  to  match  the  jet  deflection  of  the  output  stage  of 
the  preamplifier  as  closely  as  possible. 

TABLE  6.  LPA  PREAMPLIFIER  GEOMETRIC  PROPERTIES;  STAGES  1 ,  2,  AND  3 


Parameter 

Numerical  Value 

Supply  nozzle  width,  b$ 

2.5  x  10~4  m 

Supply  nozzle  height,  hg 

3.75  x  10"4  m 

Aspect  ratio,  a 

1.5 

Average  control  channel  _ 
width  normalized  by  bg,  Bc 

2.75 

Control  channel  length 
normalized  by  bs>  X£ 

10 

Discharge  coefficient,  Cd 

0.7 

Momentum  flux  discharge 
coefficient,  C 

0.65 

Blocked  load  gain,  GpB 

9.0 

Adopting  the  convention  that  stages  1,  2,  and  3  are  the  identical  pre¬ 
amplifier  stages  and  stage  4  is  the  output  amplifier  stage  makes  for  simpler 
reference  when  discussing  static  gain  and  frequency  response.  Equations  (115) 
and  (116)  are  combined  so  that  the  output  amplifier  stage,  stage  4,  can  be  matched 
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with  stage  3,  the  last  preamplifier  stage.  The  modified  Reynolds  number*  N^,  is 
held  constant,  at  its  highest  value  which  ■"'•oduces  laminar  flow,  for  all  stages 
of  the  amplification.  This  corresponds  to  a  modified  Reynolds  number  of  N^  =  120. 
Therefore,  solving  for  the  fourth  stage  aspect  ratio  results  in 

b,(3)  -y- 

3  bfw  y  y4l  *  ^27) 

The  preamplifier  parameters  are  b$(3)  3  2.5  x  10“^  m  and  o(3)  3  1.5  from 
Table  6.  The  output  amplifier  is  block  loaded  for  measurement  purposes  and  there¬ 
fore  Gp(4)  3  10.  A  standard  LPA  amplifier  with  bg(4)  3  5.0  x  10‘4  m  is  employed 
in  the  output  amplifier,  and  the  aspect  ratio  of  this  stage  is  calculated  to  be 
o(4)  3  0.24.  Viscous  losses  reduce  gain  for  amplifiers  with  an  aspect  ratio  below 
0.3;  therefore,  some  dynamic  range  is  sacrificed  to  make  o(4)  >  0.3.  The  standard 
thickness  of  the  laminations  require  that  the  minimum  aspect  ratio  of  the  fourth 
stage  be  o(4)  =  0.5.  If  the  fourth  stage  is  used  to  drive  a  load,  G  (4)  would 

r 

decrease  because  of  the  reduced  output  resistance,  and  the  optimum  value  of  cr{4) 
would  have  to  increase.  This  situation  results  in  a  better  match  between  the 
third  and  fourth  stages. 

Now  that  the  gain  block  is  defined,  the  overall  static  gain  can  be  calculated. 
Equation  (112)  gives  the  gain  of  an  amplifier  stage  as  a  function  of  its  blocked 
load  gain,  the  amplifier  output  resistance,  and  the  attached  load  resistance.  In 
a  cascaded  amplifier  gain  block  the  load  resistance  of  the  (N-l)  stage  equals 
the  input  resistance  of  the  (N)  stage.  For  the  prupose  of  measurement  the  final 
amplifier  stage  will  be  blocked  and  its  load  resistance  will  be  infinite.  This 
will  cause  the  gain  of  the  final  amplifier  stage  to  equal  its  blocked  load  gain. 

The  blocked  load  gain  of  each  stage  in  the  preamplifier  is  given  in  Table  6  as 
nine.  The  blocked  load  gain  of  the  final  amplifier  stage  is  given  to  be  ien.  This 
increase  is  due  to  an  improvement  in  the  design  of  the  final  amplifier  staye 
over  the  preamplifier  stage.  Therefore,  writing  the  gain  for  an  intermediate  stage 
of  amplification  within  the  gain  block  results  in 


,  SpB("-l> 

Gp  N‘^  =  rjn-i; 


(128) 


R-(N) 


The  ratio  of  the  output  resistance  to  the  input  resistance  can  be  solved,  using 
equations  (109)  and  (110). 


R0(N-1)  0.50  RS(N-1) 


■rTIN)  o’. 75  RjN)  • 

Applying  equation  (108)  to  equation  (3)  results  in  a  more  generalized 


(129) 


expression. 


R0(N-1) 
— R.j  (n) 


0.667 


2  (bs(N)  )3 
bs(h-l) 


(130) 


For  the  first  two  stages  of  amplification  within  the  preamplifier. 


R0(D 

R^W 


*  0.677  , 


(131) 


R0(2) 

RTT3T=  0-677  . 


(132) 


Stage  three  of  the  preamplifier  is  connected  to  the  final  stage  amplifier  and 
results  in  a  resistance  ratio  of 


Ro<3>  „  rTT  a(4)]2  bs<4!  3 

KjOT  -  °-677  [TtfrJ  [et(3tJ 

with  o(4)  *  0.5  and  b  (4)  =  5.0  x  10"4  m,  equation  (133)  is  solved. 


(133) 


Ro(3) 

S7W-0-59  • 

The  total  system  gain  with  the  last  stage  blocked  is  given  by 


(134) 


Gpsd)  Qpg(2)  ^pg(3)  Gpg(4) 
(1  +  0.677)2  (1  +  0.659) 


(135) 


>  /  _*  V.V.V  W  > 
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To  calculate  the  output  pressure  difference  of  the  gain  block,  one  must 
determine  the  input  pressure  difference.  The  static  output  pressure  difference 
is  therefore  the  static  input  pressure  difference  mutli plied  by  the  total  LPA 
system  gain  given  by  equation  (136).  The  input  pressure  difference  to  the  LPA 
gain  block,  for  a  displacement  of  the  end  of  the  bimorph  is  given  by  equation 
(119).  To  calculate  the  final  system  output  pressure  difference  several  inter¬ 
mediate  calculations  will  have  to  be  made  to  develop  values  for  the  A.  coefficients 

J 

which  are  used  by  equation  (119). 

The  first  stage  amplifier  is  operated  at  a  modified  Reynolds  number  of  NR  =  90 
to  insure  laminar  flow.  This  results  in  a  supply  pressure  and  a  supply  flow  given 
by  equations  (106)  and  (107)  respectively,  to  that  stage  of  amplification.  This 
is  given  by 

.2 


P$(l)  =  531  N/m 


Qc(l )  «  1.95  x  10 


■6  m" 


(137) 

(138) 


.-v 

.N 


.V 

a 


Equations  (55)  and  (56)  are  then  used  to  calculate  the  mass  entrained  flow. 

i-8 


We  =  5.21  x  10  °  kg/s 


(139) 


The  Reynolds  number  is  found  using  equation  (105).  Then  equation  (59)  is  solved 
for  the  control  channel  resistance  using  the  Reynolds  number  and  the  values  listed 
in  Table  5. 

Nr  =  500  ,  (140) 


Rc  =  6.05  x  107  ^ 
snr 


(141) 


The  vent  resistance  is  calculated  from  equation  (61)  to  be 


,v 

V 


R  =  8.16  x  108 

snr 


(142) 


v*  V  V  •  \  • 
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With  these  values  and  with  the  density  of  air  p,  equal  to  1.2  kg/m  ,  the  nozzle 

-4  -3 

radius,  r,  equal  to  4.32  x  10  m,  the  nozzle  length,  Lb,  equal  to  6.35  x  10  m, 

13 

the  nozzle  entrance  coefficient,  Cg  =  0.85  ,  and  the  nozzle  underlap,  u  =  3.81  x 
-5 

10  m,  one  can  calculate  the  A.  coefficients.  These  coefficients  are  solved 

J 

using  equations  (71)  through  (76)  and  are  given  by 


1.92  x  10"3  /*SL 
i  m 


A2  *  1.47  x  10 


-9 


m-s 


A,  *  5.21  x  10‘8  & 
3  s 

A.  *  5.86  x  1<T7  m2 


A5  -  1.13  x  10 
Ag  -  2.31  x  10 


-5 


/kgwri 


-3 


m 


(143) 

(144) 

(145) 
046) 

(147) 

(148) 


The  supply  pressure  to  the  bimorph,  Psp  must  now  be  calculated  to  deliver 
an  input  to  the  gain  clock  which  will  maximize  the  system  sensitivity  but  will 
not  saturate  the  LPA  gain  block.  To  do  this  a  bimorph  supply  pressure,  Pgp  is 
assumed.  Equation  (81)  is  used  to  derive  the  first  stage  amplifier  bias  pressure 
for  the  assumed  bimorph  supply  pressure.  The  intermediate  value  of  P  is  then 
calculated  with  equation  (91).  The  expression  relating  the  pressure  difference 
at  the  nozzles  to  the  deflection  of  the  end  of  the  bimorph  is  then  given  by 
equation  (119).  The  steady-state  deflection  of  the  end  of  the  bimorph  for  an 
applied  voltage  is  then  given  by  solving  equation  (97)  in  the  steady-state. 

This  results  in 


X  = 


‘VA  C5 
1  *  C0C4 


(149) 
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B.W.  Anderson,  The  Analysis  and  Design  of  Pneumatic  Systems,  John  Wiley  and 
Sons,  Inc.,  New  York,  (1976). 


The  values  of  the  Ci  coefficients  can  be  determined  using  equations  (98)  through 
(103).  Substituting  for  VA  the  maximum  voltage  to  be  delivered  to  the  bimorph 
results  in  its  maximum  displacement.  Using  equation  (119)  one  can  determine  the 
maximum  pressure  output  of  the  nozzles.  Multiplying  this  by  the  LPA  system  gain 
results  in  an  expression  for  the  maximum  output  pressure  difference  of  the  system 
for  the  maximum  input  voltage.  This  pressure  difference  is  then  checked  against 
equation  (120)  in  an  optimally  designed  gain  block  to  determine  if  the  gain  block 
has  saturated.  If  the  gain  block  is  not  optimally  designed,  the  possibility  exists 
that  an  intermediate  stage  could  saturate  prematurely.  If  the  gain  block  has 
saturated  or  is  significantly  below  saturation,  the  bimorph  supply  pressure  is 
adjusted  up  or  down  to  compensate  for  this.  This  iterative  process  is  repeated 
until  a  satisfactory  supply  pressure  to  the  bimorph  is  determined. 

This  technique  is  applied  to  the  LPA  gain  block  already  developed.  After 
several  iterations,  the  supply  pressure  to  the  bimorph  is  calculated  to  be 

Psl  *  36.0  N/m2  (2.7  x  10”1  mm  Hg)  .  (150) 

This  results  in  a  first  stage  bias  pressure,  calculated  with  equation  (81)  of 

PQ  =  34.0  N/m2  (2.56  x  10"1  mm  Hg)  .  (151) 

The  intermediate  value  of  P  is  then  calculated  using  equation  (91) 

P  =  1.41  /  N/m2  (152) 

The  expression  relating  the  pressure  difference  output  of  the  nozzles  for  an  input 
motion  of  the  end  of  the  bimorph  is  given  by  equation  (119)  and  is  shown  below. 


=  2.0  x  105 

s*  nr 


(153) 


The  maximum  steady-state  deflection  of  the  end  of  the  bimorph  for  a  maximum 
steady-state  applied  voltage  of  14  V  is  given  by  equation  (149).  To  solve  this 
equation  the  values  of  Table  5  are  substituted  into  equations  (95),  (96)  and  (98) 
throuqh  (103)  to  solve  for  the  C.  frtoffirionfe  wifh  O  «•  fi 
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C„  .  1.06  x  10-'5  . 
C2  =  3.17  x  10'5^  , 


C4  *  1.20  x  10"3  g  , 
C5  *  -3.36  x  IQ"7  m/v 
Cfi  *  0.0  s3 
C7  *  2.28  x  10“8  s2  , 
Cg  »  6.01  x  10‘6  s  . 


(154) 

(155) 

(156) 

(157) 

(158) 

(159) 

(160) 
(161) 


Solving  equation  (149)  results  in 

x|  *  5.0  x  10“6  m  .  (162) 

•max 

Substituting  equation  (162)  into  equation  (153)  results  in  an  expression  for  the 
maximum  pressure  difference  output  of  the  nozzles  for  the  maximum  input  voltage  of 
15  V. 


(Pi 


P2> 


max 


With  a  total  LPA  system  gain  of  Gp^ 
the  maximum  output  voltage  is  given 


1.01  N/m2  (7.6  x  10~3  mm  Hg) 

*  1650,  the  output  of  the  LPA  gain  block  for 
by 


(p10(4)  -  P20(4)l  *  1.66  x  103N/m2  (12.5  mm  Hg)  .  (163) 

With  the  fourth-stage  supply  pressure  at  Pg(4)  =  2.13  x  103  N/m2  (16  mm  Hg) 
calculated  from  equation  (106),  this  much  pressure  recovery  should  cause  the 
final-stage  amplifier  to  operate  in  its  nonlinear  input-output  region.  However, 
the  experimental  results  which  follow  show  that  this  is  not  the  case.  For  this 
particular  LPA  stage  there  is  slightly  better  than  70%  pressure  recovery,  but 
this  should  not  be  assumed  in  future  designs. 
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Now  that  the  system  has  been  designed  statically,  it  should  be  checked 
dynamically.  Equation  (122)  is  used  to  check  the  frequency  response  of  the 
amplifiers.  For  the  preamplifier,  with  Z  =  3,  the  upper  limit  of  the  usable 
frequency  range  of  the  final  stage  amplifier,  Z  =  1,  is  found  to  be  2400  Hz. 

The  natural  frequency  of  the  bimorph  is  given  by  equation  (14)  and  is  calculated 
to  be  1050  Hz.  To  determine  other  frequencies  that  might  give  distortions  in  the 
mgnitude  response,  equation  (136)  is  used.  Applying  this  equation  to  the  final 
stage  amplifier  shows  that  a  resonant  peak  might  be  expected  at  *  620  Hz  be¬ 
cause  of  internal  acoustic  feedback.  This  resonant  peak  is  later  noted  in  the 
experimentally  determined  system  frequency  response.  The  magnitude  of  the  peak 
is  unknown,  and  this  is  an  aspect  of  the  LPA  which  is  now  under  additional  in¬ 
vestigation^’^®. 

The  results  of  this  dynamic  check  suggest  that  a  reasonably  flat  magnitude 
response  out  to  just  below  the  resonant  frequency  of  the  bimorph  can  be  expected 
from  this  system.  The  external  series  resistor,  RA>  can  be  changed  to  vary  the 
system's  frequency  response.  This  may  help  to  further  extend  the  usable  frequency 
range  of  the  valve  system. 

A  final  point  to  be  considered  in  the  analytical  design  is  if  the  LPA  gain 
block  is  necessary  or  if  it  is  possible  to  increase  the  bimorph  chamber  supply 
pressure  until  the  same  pressure  difference  is  available  downstream  of  the  nozzles, 
as  was  available  at  the  output  of  the  LPA.  The  gain  block  is  necessary  because 
the  bimorph  is  very  limited  in  its  output  force  capability.  Any  significant 
bimorph  chamber  supply  pressure  will  cause  the  electromechanical  system  to 
exhibit  a  poorer  frequency  response  because  of  the  destabilizing  flow  forces. 

17T.M.  Drzewiecki,  A  Fluidic  Voice  Communication  System  and  Data  Link,  Ph.D. 

Thesis,  Naval  Post  Graduate  School,  Monterey,  Ca.  (March  1980). 
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T.M.  Drzewiecki,  A  High-Order,  Lumped-Parameter  Jet-Dynamic  Model  for  the 
Frequency  Response  of  Laminar  Proportional  Amplifiers,  20th  Anniversary  of 
Fluidics  Symposium,  The  American  Society  of  Mechanical  Engineers  (1980). 
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This  can  be  shown  analytically  if  the  characteristic  equation  of  equation  (97)  is 


examined.  This  is  given  by 
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Equation  (164)  is  nonlinear  with  respect  to  the  bimorph  supply  pressure,  P 
A  root  locus  can  be  plotted  for  equation  (164)  by  varying  P$^,  calculating  all 
the  intermediate  steps  and  finding  the  roots  of  equation  (164).  This  root  locus 
is  shown  in  Figure  17  for  R  =  0.0.  Figure  16  shows  that  as  P$^  is  increased, 
there  is  a  destabilizing  effect  on  the  system's  response.  For  this  reason,  the 
bimorph  supply  pressure  is  chosen  to  be  P$^  *  36.0  N/m  (2.7  x  10  mm  Hg),  which 
results  in  an  insignificant  change  in  the  dynamic  behavior  specified  by  equation  (97). 


A  prototype  model  of  the  valve  system  was  then  constructed  using  the 
dimensions  outlined  in  the  analyticaldesign  section.  An  assembly  diagram  for 
the  design  is  shown  in  Figure  17.  The  prototype  was  built  attempting  to  simplify 
its  construction,  minimize  volumes  to  minimize  fluid  capacitance,  and  minimize 
passage  lengths  for  reduced  signal  propagation  time.  This  resulted  in  a  device 
which  was  inexpensive  to  build  and  which  delivered  maximum  bandwidth. 

The  assembly  diagram  of  Figure  17  shows  the1 construction  of  the  housing 
containing  the  flapper-nozzle  and  the  first  stage  of  pneumatic  amplification. 

Only  the  first  LPA  stage  is  shown  to  avoid  confusion.  The  diagram  also  contains 
lines  which  are  drawn  for  the  purpose  of  flow  visualization. 

The  two  pneumatic  inputs  to  the  flapper-nozzli  housing  are  the  supply  pressure 
to  the  bimorph  and  the  supply  pressure  to  the  LPA  gain  block.  The  flow  through 
nozzles  is  controlled  by  the  position  of  the  end  of  the  bimorph.  Figure  18  shows 
the  bimorph  deflected  with  more  flow  entering  one  nozzle  than  the  other.  The 
relative  amount  of  flow  is  indicated  by  the  thickness  of  the  line  used  to  re- 
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present  that  flow.  The  outputs  of  the  two  nozzles  are  then  channeled  up  through 
several  laminations  to  the  LPA  lamination  where  these  output  pressures  are  used 
as  the  control  input.  The  controls  deflect  the  laminar  supply  jet,  and  a  resulting 
output  pressure  difference  is  available  at  the  LPA  output  ports.  The  output  of 
tnis  stage  is  then  available  to  be  channeled  up  to  the  input  of  a  second  LPA 
staoe.  This  process  is  repeated  for  four  stages  of  amplification  in  the  working 
prototype. 

5.3  Testing  and  Verification 

The  valve  system  is  tested  both  statically  and  dynamically  to  determine 

the  accuracy  of  the  design  equations  used.  The  static  input-output  relationship 

for  the  system  is  shown  in  Figure  18.  The  system  has  a  linear  operating  region 

?  2 

of  +  7.0  volts  about  0.0  volts  with  a  system  sensitivity  of  1.2  x  10  N/m  -V, 

(0.9  mm  Hg/V).  This  result  compares  favorably  with  the  expected  sensitivity  of 
1.1  x  102  N/m2-V,  (0.8  mm  Hg/V).  The  expected  sensitivity  is  derived  by 
dividing  the  maximum  output  pressure  difference  of  the  gain  block  given  by 
equation  (163)  by  the  maximum  input  voltage  of  +  7.0  volts. 

Figure  18  also  demonstrates  that  the  system  output  exhibits  hysteresis. 

The  hysteresis  is  an  inherent  characteristic  of  the  bimorph  caused  by  the  polari¬ 
zation  of  the  crystal  resulting  from  the  applied  supply  voltage.  The  cause  of 
this  effect  is  explained  earlier  in  Section  2. 

The  frequency  response  of  the  system  was  measured  for  a  bender  series  re¬ 
sistance  of  Ra  *  0.0  ohm.  The  results  of  that  test  along  with  the  analytically 
predicted  response  of  the  system  are  shown  in  Figure  19.  There  are  some  incon¬ 
sistencies  between  the  measured  and  the  predicted  response  on  the  magnitude  plot. 
The  most  noticeable  inconsistency  is  the  difference  in  the  size  of  the  resonant 
peaks.  The  measured  resonant  peaks  are  smaller  because  the  pneumatic  input  to 
the  preamplifier  is  so  large,  for  the  electrical  input  magnitude  chosen,  that 


Figure  18.  Static  Input-Output  Relationship  of 
Valve  System 

1  Large  Division  Horizontal  -  2.9V 

2 

1  Large  Division  Vertical  *  410  N/m 


Figure  19.  Measured  and  Predicted  System  Frequency  Response 
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the  preamplifier  saturates.  This  saturation  occurs,  by  design,  for  relatively 

small  electrical  inputs  at  the  resonant  frequency  to  obtain  the  maximum  gain  and 

bandwidth  in  the  operating  frequencies.  A  way  to  reduce  this  effect  would  be 

to  place  the  bender  resonance  beyond  the  amplifier  bandwidth  limits.  The  LPA 

would  then  filter  the  resonant  frequency  of  the  bender  and  eliminate  the  saturation. 

Figure  19  demonstrates  that  the  system  amplitude  versus  freqeuncy  response 

exhibits  a  resonance  at  650  Hz.  This  resonance  is  a  result  of  internal  acoustic 

feedback  in  the  fourth-stage  amplifier.  To  analytically  predict  the  frequency  of 

maximum  gain  due  to  this  phenomenon,  equation  (126)  is  utilized.  With  the  amplifier 

-3 

nozzle  to  splitter  distance,  Xsp,  equal  to  4.0  x  10  m,  the  nozzle  discharge  co- 

3 

efficient,  C^,  equal  to  0.7,  the  amplifier  supply  pressure,  Ps»  equal  to  4.8  x  10 

3 

PA,  and  the  density  of  air,  p,  equal  to  1.2  kg/m  ,  the  frequency  of  the  jet 
oscillation,  1/T,  is  calculated  with  equation  (126)  to  be  620  Hz.  This  compares 
favorably  with  the  650  Hz  resonance  peak  of  Figure  19. 

A  comparison  of  the  phase  angle  versus  frequency  plots  of  Figure  20  shows  good 
agreement  between  the  measured  response  and  the  calculated  response.  If  90° 
phase  lag  is  the  cutoff  point  for  the  usable  frequency  range  of  the  system,  this 
system  has  a  usable  frequency  range  from  a  DC  level  to  140  Hz. 

If  magnitude  information  is  important  in  a  system  design,  then  sending  a 

constant  amplitude  signal  within  the  bandwidth  of  this  device  will  result  in 
minimum  magnitude  distortion  of  the  input  signal.  This  results  because  the 
system  phase  lag  is  due  almost  entirely  to  a  transport  delay  due  to  the  pneumatic 

signal  path  length.  A  series  resistor  may  be  added  to  the  bender  electrical  supply 

to  flatten  the  magnitude  response  over  a  greater  frequency  range.  Figure  20. 
reveals  the  measured  and  predicted  response  for  a  series  resistor,  =  0.2  inn. 

The  plot  shows  that  for  a  tolerance  of  _+  3  db  of  magnitude  distortion,  the  system 
exhibits  a  nearly  flat  amplitude  versus  frequency  response  and  a  linear  phase 
angle  versus  frequency  relationship  up  to  1000  Hz. 
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The  test  data  indicate  that  the  valve  design  does  work.  The  data  also 
indicate  that  the  static  and  dynamic  results  are  close  to  those  predicted  by  the 
design  equation.  Now  that  a  working  prototype  has  been  built,  obvious  improvements 
can  be  made.  One  improvement  is  to  reduce  the  pneumatic  signal  path  length. 

This  should  shorten  the  pneumatic  time  delay  and  improve  the  system's  phase 
response.  Another  improvement  would  be  to  reduce  the  dimensions  of  the  bimorph. 

This  would  improve  the  magnitude  response  of  the  system.  A  shorter  bimorph  would 
result  in  a  reduced  motion  at  the  end  of  the  bimorph  for  an  applied  voltage.  This 
can  be  compensated  for  with  an  improved  amplifier  system.  The  amplifier  system 
can  be  designed  so  that  all  stages  operate  at  a  high  supply  pressure.  A  high  supply 
pressure  results  in  a  high  supply  velocity.  Equation  (122)  shows  that  a  gain  block 
with  a  high  supply  velocity  will  exhibit  improved  magnitude  versus  frequency 
response.  A  high  supply  velocity  also  will  reduce  the  time  required  for  pneumatic 
signal  propagation  through  the  gain  block.  This  should  further  improve  the  phase 
response  of  the  system.  The  final  section  concerning  future  considerations, 
analaytically  develops  a  design  that  snould  exhibit  an  improved  response. 

6.  FUTURE  CONSIDERATIONS 

Much  can  be  learned  about  a  design  once  an  initial  prototype  has  been  built. 

Now  that  the  original  design  has  been  shown  to  work  and  the  design  equations 
correctly  predict  its  response,  a  new  design  can  be  developed  using  those 
equations.  The  new  design  should  extend  the  usable  frequency  range  of  the  device 
through  a  better  choice  of  the  piezoelectric  bimorph  and  the  attached  LPA  gain 
block.  A  summary  of  the  important  physical  properties  of  the  new  prototype  valve 
valve  system  is  provided  in  Appendix  B. 

The  new  design  objective  is  to  improve  the  frequency  response  of  the  system. 
This  is  accomplished  by  reducing  the  dimensions  of  the  bimorph  and  the  LPA. 

Again,  the  width  and  thickness  of  the  bimorph  is  fixed  due  to  the  specified 
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minimum  nozzle  diameter  (0  *  8.6  x  10  m)  and  the  availability  of  standard  bimorph 
sizes.  Therefore  the  only  remaining  bimorph  dimension  that  may  be  picked  is  its 
length.  The  length  of  the  bimorph  should  be  reduced  as  much  as  possible  to  maximize 
the  frequency  response.  However,  too  short  a  bender  will  cause  the  nozzle  output 
pressure  signal  to  be  so  small  that  detecting  it  with  the  first  stage  amplifier 
will  be  difficult.  To  calculate  a  reduced  length  for  the  bimorph,  a  first-stage 
amplifier  must  be  assumed.  Then,  with  a  knowledge  of  the  gain  and  dynamic  range 
of  the  first  stage  amplifier,  one  can  calculate  a  suitable  smaller  bimorph  length. 
This  situation  requires  a  certain  feel  for  what  might  work  best  in  the  first- 
stage  amplifier. 

Because  it  is  desired  to  improve  the  frequency  response  of  the  system,  the 

-4 

first  stage  amplifier  is  chosen  to  have  a  supply  nozzle  width  of  b$(l)  =  2.5  x  10  m 
and  an  aspect  ratio  of  o(l)  *  1.0.  This  results  in  a  maximum  supply  pressure  of 
Ps(l)  *  2.13  x  103  N/m2,  (16  mm  Hg)  from  equation  (106).  The  usable  frequency 
range  of  this  amplifier  has  been  increased  because  of  the  choice  of  the  aspect 
ratio  and  the  supply  nozzle  width.  This  maximum  usable  frequency  can  be  calculated 
with  equation  (122)  and  is  given  by  f  =  4800  Hz. 

To  determine  the  smallest  length  of  a  bimorph,  which  would  produce  a  de¬ 
tectable  input  pressure  signal  to  the  first-stane  amplifier,  one  notes  from 
equation  (22)  that  the  deflection  of  the  end  of  the  bimorph  varies  with  the  square  of 
its  length.  Using  the  same  valve  underlap  in  the  design  described  here  as  was 
used  in  the  first  design  found  in  Section  5  (u  ®  3.8  x  10"^  m),  the  pressure 
output  of  the  valve  varies  proportionally  with  the  deflection  of  the  bimorph; 
therefore,  the  pressure  output  of  the  flapper-nozzle  configuration  must  also  vary 
with  the  square  of  the  bimorph's  length.  With  this  knowledge  of  the  first  stage 
IPA  geometry  and  the  relationship  between  the  bimorph's  length  and  the  output 
pressure  difference  available  at  the  nozzles,  a  new  bimorph  length  can  be 
calculated. 
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Equation  (119)  states  that  the  linear  pressure  recovery  for  an  LPA  stage  is 

3  o 

70  percent  of  that  stage's  supply  pressure.  With  P$(l)  *  2.13  x  10  N/m  (16  mm  Hg) 
this  corresponds  to  a  maximum  pressure  recovery  of  1.49  x  10  N/m  (11.2  mm  Hg). 

It  will  be  assumed  that  the  first  stage  operates  with  a  gain  of  10.  This  may  be 
high,  because  the  output  is  not  blocked,  but  it  will  give  a  good  approximation. 

2 

This  means  that  the  maximum  input  pressure  difference  to  this  stage  is  1.49  x  10 
N/m  (1.12  mm  Hg).  Drzewiecki  has  shown  that  a  dynamic  range  of  100,000:1  can  be 
expected  for  an  LPA  stage  of  this  design^.  Since  the  dynamic  range  is  a  ratio  of 
the  largest  input  to  the  smallest  input  which  gives  uniform  gain,  the  smallest 
allowable  input  pressure  difference  to  the  first-stage  amplifier  can  now  be 
calculated.  This  is  given  by  1.5  x  10  N/m  (1.1  x  10  mm  Hg).  Because  some 
dynamic  range  in  the  electrical  input  signal  is  desired,  the  smallest  allowable 

I 

input  pressure  difference  to  the  first  stage  is  arbitrarily  increased  to  5.0  x 
10" 1  N/m2  (3.8  x  10"3  mm  Hg)/ 

Equation  (163)  gives  an  expression  for  the  maximum  output  pressure 
difference  available  at  the  nozzles,  for  a  maximum  input  voltage  of  15  V.  This 

expression  is  for  the  working  prototype  of  Section  5  with  a  corresponding  length, 
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L  *  1.5  x  10  m,  and  a  corresponding  output  pressure  difference  of  1.01  N/m 

-3 

(7.6  x  10  mm  Hg).  Since  the  square  of  the  length  of  the  proposed  prototype 
should  vary  with  the  pressure  output  available  downstream  of  the  nozzles,  the 
characteristics  of  the  working  prototype  can  be  used  to  estimate  a  minimum  length 
of  the  bimorph  for  the  proposed  prototype.  This  corresponds  to  L  -  1.05  x  10"^  m. 
The  new  piezoelectric  bimorph  physical  and  material  properties  are  summarized  in 
Table  7. 

The  equations  of  Section  3  can  then  be  used  to  calculate  the  electromechanical 
system  parameters  in  the  same  manner  that  they  were  calcualted  in  that  section. 

^T.M.  Drzewiecki,  A  Fluidic  Voice  Communication  System  and  Data  Link,  Ph.D. 

Thesis,  Naval  Post  Graduate  School,  Monterey,  Ca.  (March. 1980). 


TABLE  7.  PHYSICAL  AND  MATERIAL  PROPERTIES  OF  LEAD  ZIRCONATE  TITANATE 
CERAMIC  BIMORPH  USED  FOR  PROPOSED  PROTOTYPE 


Property 

Measured  Value 

Length,  L 

1 .05  x  10‘2  m 

Width,  W 

1.4  x  10"3  m 

Thickness,  T 

4.9  x  10-4  m 

Nickel  plate  thickness,  t^ 

1.52  x  10*6  m 

Lead  zirconate  titanate  thickness,  t2 

1.93  x  10"4  m 

Brass  shim  thickness,  t3 

1.03  x  10”4  m 

Deni  sty,  p 

7.5  x  103  kg/m3 

Piezoelectric  charge  coefficient,  d.^ 

-1.8  x  10-10  m/V 

Piezoelectric  voltage  coefficient,  g^-j 

-1.1  x  10'2  V-m/N 

Nickel  elastic  modulus,  E1 

1.26  x  1011  N/m2 

Lead  zirconate  titanate 
elastic  modulus,  E2 

1.09  x  1011  N/m2 

Brass  elastic  modulus,  E^ 

9.0  x  1010  N/m2 

Those  results  are  summarized  in  Table  8.  It  can  be  noted  from  these  parameters 
that  the  natural  frequency  of  the  bimorph  has  increased  to  2300  Hz.  These  results 
are  useful  later  in  calculating  an  improved  flapper-nozzle  supply  pressure. 

Since  the  first  stage  of  the  gain  block  has  already  been  designed  and  an 
approximate  input  sensitivity  to  that  stage  has  been  determined,  it  is  logical  to 
finish  the  design  of  the  gain  block.  With  that  completed,  the  supply  pressure 
to  the  flapper  nozzle,  which  gives  a  maximum  output  pressure  difference,  can  be 
calcualted  in  the  same  manner  it  was  done  in  Section  5. 

The  final  three  stages  of  amplification  are  chosen  to  be  self-staged.  Even 
though  self-staging  an  amplifier  three  times  results  in  a  reduced  usable  frequency 


TABLE  8.  ANALYTICALLY  PREDICTED  SYSTEM  PARAMETERS  USED  FOR 
PROPOSED  PROTOTYPE 


Parameter 

Analytically  Predicted  Value 

System  spring  constant,  K$ 

2636  N/m 

Effective  mass,  M 

1.27  x  10"5  kg 

System  damping,  B 

5  x  10‘3  N-s/m  (estimated) 

Piezoelectric  force  constant,  C^ 

-4.11  x  10"4  N/V 

Piezoelectric  back  current 
constant, 

-3.43  x  10'4  A-s/m 

Bimorph  capacitance,  Cj 

5.0  x  10"10  F 

range,  the  frequency  range  of  each  stage  is  so  high  that  the  reduced  range  is 
still  very  great.  The  aspect  ratio  of  stages  two,  three,  and  four  are  given  to  be 
o( 2)  *  a( 3)  *  a(4)  »  0.6.  The  supply  nozzle  width  for  each  stage  is  given  to  be 
bs(2)  =  b$(3)  =  b$(4)  =  2.5  x  10*4  m.  Equation  (106)  defines  the  supply  pressure 
to  these  stages  as  P$ ( 2)  *  P$(3)  =  Ps(4)  =  5.8  x  103  N/m2  (44  mm  Hg).  This 
amplifier  design  gives  improved  magnitude  response,  with  Z  =  3,  as  f  *  4000  Hz. 

The  total  gain  of  this  LPA  gain  block  can  be  calcualted  using  the  LPA  design 
equations  of  Section  4  in  the  same  manner  as  they  were  used  to  calculate  the 
total  gain  of  the  gain  block  in  Section  5.  This  results  in  a  total  system  gain, 
with  the  last  stage  blocked,  given  by 

Gpt  =2900  (165) 

Multiplying  the  total  gain  of  the  system,  given  in  equation  (165)  by  the 

available  input  pressure  difference  to  the  first  stage  given  by  5.0  x  10*^  N/m2 
-3  . 

(3.8  x  10  mm  Hg),  results  in  a  maximum  output  pressure  difference  of  the  LPA 
gain  block  given  by  1.16  x  103  N/m2  (8.7  mm  Hg).  This  value  is  well  within  the 
linear  pressure  recovery  of  the  last  stage  amplifier,  and  the  design  should  work. 
To  calculate  a  better  value  for  the  flapper-nozzle  supply  pressure,  the 


iterative  process  of  Section  5  is  used.  After  several  iterations,  the  flapper- 
nozzle  supply  pressure,  which  gives  a  maximum  of  70-percent  pressure  recovery 
to  any  intermediate  stage  within  the  gain  block,  is  calculated  to  be 

P$1  =  290  N/m2  (2.2  mm  Hg)  (166) 

This  results  in  a  maximum  steady-state  output  pressure  difference  of  the  LPA 
gain  block  for  the  maximum  input  voltage  of  15  V  given  by 

pQ1(4)  -  PQ2(4)j  =  4-°  x  1q3  N/m2  (3-°  x  1q1  h9)  (167) 

•max 

The  increased  flapper-nozzle  supply  pressure  will  not  affect  the  stability  of  the 
system.  This  can  be  determined  by  evaluating  the  roots  of  the  system's 
characteristic  equation  given  by  equation  (164)  and  noting  that  they  are  in  the 
left  half  plane. 

The  design  presented  in  this  section  should  exhibit  an  improved  frequency 
response  and  greater  output  sensitivity  over  the  design  presented  and  tested  in 
Section  5.  The  improved  frequency  response  is  due  to  a  shorter  bimorph  length  and 
an  LPA  gain  block  engineered  for  fast  response.  The  increased  gain  is  a  result 
of  a  better  impedance  match  between  individual  LPA  stages  within  the  gain  block. 
The  final  result  is  an  improved  system  response. 

A  final  consideration  for  a  complete  servovalve  design  is  a  last  stage  power 
amplifier.  The  power  output  of  the  LPA  is  somewhat  limited  by  its  small  size 
and  the  requirement  that  it  run  with  a  relatively  low  supply  pressure  so  that  the 
jet  remains  laminar.  Any  type  of  power  amplifier  that  mignt  be  added  would  have 
a  tendency  to  reduce  the  bandwidth  of  this  system.  This  will  be  true  because  a 
moving  part  power  amplifier  will  have  mechanical  parts.  This  will  require 
acceleration  forces  and  some  fluid  to  mechanical  force  amplification,  which  will 
introduce  fluid  capacitance.  An  example  of  this  would  be  some  sort  of  a  diaphragm 
amplifier.  A  relatively  large  diaphragm  may  be  necessary  for  sufficient  pneumatic 


amplification.  The  motion  of  the  diaphragm  will  introduce  a  fluid  capacitance  to 
the  system.  The  overall  result  will  be  a  system  which  is  capable  of  a  far  slower 
response  but  has  a  higher  output  power  capabilities.  An  interesting  field  of 
study,  for  future  research,  would  be  the  development  of  a  fluidic  power  amplifier 
that  would  match  an  LPA  stage.  A  power  amplifier  of  this  type  would  not  only 
have  application  in  this  project  but  also  in  the  other  areas  of  LPA  research 
where  the  sensing  capabilities  of  the  LPA  are  used  in  a  computer  control  scheme. 

|Q 

This  includes  such  projects  as  temperature  snesing  and  control  or  the  use  of 

20 

the  laminar  jet  angular  rate  sensor  in  conjunction  with  the  LPA  in  a  complete 
computer  controlled  closed  loop  system. 
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Cr  ~  1/ 


rtc, 


ks7r^T 


m/V 


CtMR 

^r 

ctbr+m 

"V 

(CTKs+ClKq)R+B 

T 

Discharge  coefficient,  dimensionless 
Entrance  coefficient,  dimensionless 
Piezoelectric  capitance,  F 
LPA  supply  nozzle  momentum  flux  discharge,  dimensionless 
Nozzle  diameter,  m 
Piezoelectric  d  coefficient,  m/V 

2 

Equivalent  modulus  of  elasticity,  N/m 

2 

Nickel  modulus  of  elasticity,  N/m 

2 

Lead  zirconate  titanate  modulus  of  elasticity,  N/m 

2 

Brass  modulus  of  elasticity,  N/m 

Force  on  the  control  volume  exerted  by  the  bimorph,  N 

Maximum  usable  LPA  frequency,  Hz 

Acoustic  feedback  frequency,  s"^ 

Force  on  the  bimorph  exerted  by  the  fluid,  N 

Force  to  deflect  bimorph  a  distance  a,  N 

Piezoelectric  g  coefficient,  V-m/N 

LPA  gain,  dimensionless 

LPA  blocked  load  gain,  dimensionless 

Total  LPA  system  gain,  dimensionless 

LPA  supply  nozzle  height,  m 

4 

Moment  of  area,  m 


Electrical  signal  field  direction,  dimensionless 

Integer,  dimensionless 

Stress  direction,  dimensionless 

Coupling  coefficient,  dimensionless 

Pizeoelectric  back  current  constant,  A-s/m 

System  spring  constant,  N/m 

Bimorph  length,  m 

Nozzle  length,  m 

Control  volume  length,  side  1,  m 
Control  volume  length,  side  2,  m 
Effective  mass  of  bimorph, 

Integer,  dimensionless 
Integer,  dimensionless 
Reynolds  number,  dimensionless 
Modified  Reynolds  number,  dimensionless 
Unit  normal  vector,  dimensionless 
Integer,  dimensionless 

Coefficient  used  to  simplify  notation,  /  N/i m 
Bimorph  polarity,  dimensionless 

2 

LPA  operating  point  bias  pressure,  N/m 

2 

Nozzle  1  pressure,  N/m 

2 

Nozzle  2  pressure,  N/m 

2 

LPA  operating  point  bias  pressure,  side  1,  N/m 

2 

LPA  operating  point  bias  pressure,  side  2,  N/m 
Atmospheric  pressure,  N/m 

2 

LPA  input  pressure,  side  1,  N/m 

2 

LPA  input  pressure,  side  2,  N/m 

2 

LPA  output  pressure,  side  1,  N/m 


2 

LPA  output  pressure,  side  2,  N/m 

LPA  control  channel  pressure,  side  1,  N/m 

LPA  control  channel  pressure,  side  2,  N/m 

2 

LPA  supply  pressure,  N/m 

2 

Flapper-nozzle  supply  pressure,  N/m 
2 

Vent  pressure,  N/m 
Charge  produced,  coulomb  (C) 

3 

Volumetric  jet  entrainment  flo  <,  ft  /s 

3 

Volumetric  LPA  supply  flow,  ft  /s 
Equivalent  electrical  resistance,  ohm  {a) 
Nozzle  radius,  m 
External  electrical  resistance, 

4 

LPA  control  channel  resistance,  kg-m  /s 

4 

LPA  input  resistance,  kg-m  /s 

4 

LPA  load  resistance,  kg-m  /s 

4 

LPA  output  resistance,  kg-m  /s 

4 

LPA  supply  nozzle  resistance,  kg-m  /s 
Bimorph  leakage  resistance,  a 

4 

LPA  vent  resistance,  kg-m  /s 
Laplace  variable,  s-1 
Bimorph  thickness,  m 
Time,  s 

Laminar  jet  transport  time,  s 

Nickel  plate  thickness,  m 

Lead  zirconate  titanate  thickness,  m 

Brass  shim  thickness,  m 

Nozzle  underlap,  m 

LPA  stage 


T 


U2 

3 

LPA  stage 

U3 

3 

LPA  stage 

U4 

3 

Amplifier  average  jet  particle  velocity,  m/s 

US 

= 

Amplifier  supply  velocity,  m/s 

V 

= 

Node  voltage,  V 

V1 

= 

Control  volume  fluid  velocity,  side  1,  m/s 

V2 

3 

Control  volume  fluid  velocity,  side  2,  m/s 

VA 

= 

Supply  voltage,  V 

V 

= 

Velocity  rector,  m/c 

V- 

= 

3 

Volume,  m 

W 

3 

Bimorph  width,  m 

W1 

= 

Control  volume  mass  flow  rate,  side  1,  kg/s 

Wi 

3 

Steady  state  continuity  equatin 

U2 

- 

Control  volume  mass  flow  rate,  side  2,  kg/s 

W2 

3 

Steady  state  continuity  equation 

Mli 

3 

Mass  flow  rate  into  cont, ol  volume,  side  1,  kg/s 

W2i 

3 

Mass  flow  rate  into  control  volume,  side  2,  kg/s 

W10 

3 

Mass  flow  rate  out  of  control  volume,  side  1,  kg/s 

W20 

3 

Mass  flow  rate  out  of  control  volume,  side  2,  kg/s 

WA1 

a 

Vent  flow,  side  1,  kg/s 

WA2 

3 

Vent  flow,  side  2,  kg/s 

We 

3 

Mass  jet  entrainment  flow,  kg/s 

WE1 

3 

Nickel  plate  equivalent  width,  m 

WE2 

3 

Lead  zirconate  titanate  equivalent  width,  m 

WE3 

3 

Brass  shim  equivalent  width,  m 

X 

3 

Bimorph  displacement,  m 

Xo 

3 

Bimorph  operating  point  deflection,  m 

O 

•X 

* 

Bimorph  operating  point  velocity,  m/s 

-/v 

v  .  /  /  '  7>X  -A  >f  •  / 
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X 


LPA  stage 
LPA  stage 
LPA  stage 

LPA  nozzle  to  splitter  distance,  m 
LPA  stage 

Number  of  LPA  self  staged  amplifier  stages  dimensionless 
LPA  stage 

2 

Coefficient  of  elasticity,  cm  /dyne 

2 

Coefficient  of  elasticity,  cm  /dyne 

2 

Coefficient  of  elasticity,  cm  /dyne 

Deflection  of  the  free  end  of  a  cantilever  mounted  beam  from 
an  applied  force,  m 
Absolute  dielectric  coefficient,  F/m 
Bimorph  damping  ratio,  dimensionless 
Angle,  rad 

2 

Kinematic  viscosity,  m  /s 

3 

Fluid  density,  kg/m 
LPA  aspect  ratio,  dimensionless 
Period  of  LPA  laminar  jet  oscillation,  s 
Bimorph  natural  frequency,  rad/s 
Bimorph  natural  frequency,  rad/s 
Direction  cosine,  dimensionless 
Direction  cosine,  dimensionless 
Direction  cosine,  dimensionless 


yi  V  •  yi  .r.T.VV’-’.VA 
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APPENDIX  A  —  SUMMARY  OF  THE  PROPERTIES  OF  THE  WORKING  PROTOTYPE 

A  summary  of  the  important  physical  properties  of  the  working  prototype  de¬ 
veloped  in  Section  5  is  provided  below  for  reference. 


Flapper-Nozzle  Housing 


Property 

Bimorph  length,  L 
Bimorph  width,  W 
Bimorph  thickness,  T 
Nozzle  underlap,  u 
Nozzle  diameter,  D 
Nozzle  length, 

Nozzle  entrance  coefficient,  Cg 
Bimorph  chamber  supply  pressure,  P^ 
Air  density,  p 


LPA  Gzin  Block 

Stages  1,  2,  and  3: 

Property 

Supply  nozzle  width,  b$ 

Supply  nozzle  height,  h$ 

Aspect  ratio,  o 

Average  control  channel  width  normalized  by 
bs*  'Bc 

Control  channel  length  normalized  by  b$,  Xc 

Discharge  coefficient, 

Momentum  flux  discharge  coefficient,  Cf 


Blocked  load  pressure  gain,  GpB 
Supply  pressure,  P.(l),  P_(2),  F . ( 3 ) 


Measured  Value 

1.5  x  10“Z  m 
1.4  x  10*3  m 
4.9  x  10'4  m 
3.8  x  10"5  m 

8.6  x  10‘4  m 
6.35  x  10'3  m 
0.85 

3.6  x  101  N/m2 
1.2  kg/m3 


Measured  Value 
2.5  x  10‘4  m 


3.75  x  10"4  m 
1.5 

2.75 

10 

0.7 

0.65 

9.0 

5.31  x  102  N/m2 


/ 


Stage  4: 

Property 

Supply  nozzle  width,  t>s 

Supply  nozzle  height,  h$ 

Aspect  ratio,  a 

Average  control  channel 
Width  normalized  by  bs,  Bc 

Control  channel  length  normalized  by  b  ,  Xc 

Discharge  coefficient,  Cd 

Momentum  flux  discharge  coefficient,  CQ 

Blocked  load  pressure  gain,  Gpe 

Supply  pressure,  P$(4) 


Measured  Value 
5.0  x  10"4  m 
2.5  x  10'4  m 
0.5 
2.75 

10 

0.7 

0.65 

9.0 

2.13  x  103  N/m2 


Figure  A-l  is  a  schematic  drawing  of  the  working  prototype.  This  figure 
illustrates  the  flapper-nozzle  housing  along  with  four  stages  of  pneumatic 
amplification.  The  maximum  static  output  pressure  difference  of  the  nozzles  is 
given  for  the  maximum  static  displacement  of  the  end  of  the  bimorph.  This 
pressure  difference  is  then  applied  across  the  input  port  of  an  LPA  gain  block. 
The  important  geometric  properties  along  with  the  cascaded  gain  and  the  resulting 
output  pressure  difference  of  each  LPA  stage  are  given. 


Figure  A-T.  Schematic  Drawing  of  Working  Prototype 


APPENDIX  B  ~  SUMMARY  OF  THE  PROPERTIES  OF  THE  PROPOSED  PROTOTYPE 


This  appendix  provides  a  summary  of  the  important  physical  properties  of 

proposed  prototype  developed  in  Section  5. 

Flapper-Nozzle  Housing 

Property 

Measured  Value 

Bimorph  length,  L 

1.05  x  10~2  m 

Bimorph  width,  W 

1.4  x  10"3  m 

Bimorph  thickness,  T 

4.9  x  10~4  m 

Nozzle  underlap,  u 

3.8  x  10"5  m 

Nozzle  diameter,  D 

•4 

8.6  x  10  m 

Nozzle  length,  Lb 

6.35  x  10“3  m 

Nozzle  entrance  coefficient, 

0.85 

Bimorph  chamber  supply  pressure,  P$-| 

2.90  x  102  N/m2 

Air  density,  p 

1.2  kg/m3 

LPA  Gain  Block 

Stages  1: 

Property 

Measured  Value 

Supply  nozzle  width,  bs 

2.5  x  10'4  m 

Supply  nozzle  height,  h$ 

2.5  x  10'4  m 

Aspect  ratio,  a 

1.0 

Average  control  channel  width  normalized  by 

b  ,  B 
s’  c 

2.75 

Discharge  coefficient, 

0.7 

Momentum  flux  discharge  coefficient,  C0 

0.65 

Blocked  load  pressure  gain,  Gpg 

10.0 

Supply  pressure,  P$(l) 

2.13  x  102  N/m2 
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Stages  2,  3,  and  4: 

Property 

Supply  nozzle  width,  bg 

Supply  nozzle  height,  hg 

Aspect  ratio,  a 

Average  control  channel  __ 

Width  normalized  by  b  ,  B 

j  L 

Control  channel  length  normalized  by  b  ,  X 

5  w 

Discharge  coefficient,  Cd 
Momentum  flux  discharge  coefficient,  CQ 
Blocked  load  pressure  gain,  Gpg 
Supply  pressure,  P$ (2) ,  Ps(3),  P$(4) 


Measured  Value 

2.5  x  10"4  m 

1.5  x  10"4  m 
0.6 

2.75 

10 

0.7 

0.65 

10.0 

5.8  x  103  N/m2 


Figure  B-l  is  a  schematic  drawing  of  the  proposed  prototype.  This  figure 
illustrates  the  flapper-nozzle  housing  along  with  four  stages  of  pneumatic 
amplification.  The  maximum  static  output  pressure  difference  of  the  nozzles  is 
given  for  the  maximum  static  displacement  of  the  end  of  the  bimorph.  This 
pressure  difference  is  then  applied  across  the  input  port  of  an  LPA  gain  block. 
The  important  geometric  properties  along  with  the  cascaded  gain  and  the  resulting 
output  pressure  difference  of  each  LPA  stage  are  given. 
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Figure  B-l.  Schematic  Drawing  of  Proposed  Prototype 
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